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ABSTRACT 


New  short-period  data  on  source  functions  both  of  explosions  and 
earthquakes  are  reported.  Studies  of  wave  propagation  in  inhomo¬ 
geneous  regions  of  the  earth  are  described.  First  results  from 
NORSAR  are  presented.  Long-period  noise  and  signal  observa¬ 
tions  have  been  made,  and  some  idea  of  NORSAR* s  discriminative 
capability  is  emerging.  Network  approaches  to  discrimination  are 
being  expanded.  Major  improvements  in  P-wave  detection  capa¬ 
bility  are  in  sight,  with  joint  array  and  network  analysis.  New 
programs  are  being  applied  to  the  rapid  acquisition  of  library  in¬ 
formation  on  seismic  events.  An  increasing  use  of  LASA  for  other 
geophysics  projects  is  apparent. 
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SUMMARY 


This  is  the  sixteenth  Semiannual  Technical  Summary  of  the  Seismic  Discrimination  Group 
of  Lincoln  Laboratory.  Our  work  is  directed  toward  acquiring  as  extensive  an  understanding 
as  possible  of  seismic  techniques  for  detecting,  locating  and  identifying  earthquakes  and  under¬ 
ground  explosions,  with  a  view  to  providing  technical  background  necessary  to  any  discussions 
on  a  treaty  banning  nuclear  weapons  tests.  To  the  present,  our  work  has  been  predominantly 
centered  on  the  most  effective  use  of  LASA,  the  large  aperture  seismic  array  in  Montana.  Our 
interests  are  now  broadening  into  the  use  of  two  recently  completed  arrays,  NORSAR  and  ALPA, 
the  World-Wide  Standard  Seismographic  Network,  and  new  long-period  instruments.  Wre  con¬ 
tinue  to  pursue  research  on  the  phenomenology  of  the  seismic  source  and  on  several  features 
of  wave  propagation.  An  increasing  amount  of  pure  geophysics  is  done  with  LASA,  and  a  por¬ 
tion  of  this  is  described  in  Sec.  VII  of  this  report. 

LASA  and  NORSAR  have  both  been  used  to  study  explosions  in  order  to  try  and  improve  our 
understanding  of  teleseismic  aspects  of  the  explosive  source.  P-waves  at  NORSAR  from  Nevadan 
explosions  for  a  wide  range  of  yields  and  burial  depths  have  been  studied  in  order  to  attempt  to 
isolate  effects  arising  from  the  primary  pulse,  the  surface  reflection  (pP),  and  the  spalling  or 
slap-down  that  occurs  shortly  afterward.  Very  clear  contributions  from  each  of  these  phenomena 
can  be  isolated.  This  study  is  corroborated  by  studies  of  core  phases  from  Cannikin  and  Milrow 
in  the  Aleutians  from  which  similar  conclusions  emerge.  It  has  further  been  possible  to  begin 
making  estimates  of  an  earthquake's  size  using  short-period  data.  Core  phases  have  been  used 
to  obtain  an  unencumbered  view  of  the  source,  and  remarkably  short  durations  for  earthquakes 
are  measured.  These  convert  into  sizes  for  events  of  less  than  10  km  even  up  to  m^6.0. 

In  recent  years,  large  computers  have  allowed  new  approaches  to  problems  in  wave  prop¬ 
agation.  One  such  approximate  technique  has  been  developed  which  handles  inhomogeneous 
regions  at  a  relatively  low  expenditure  of  computing  time.  We  report  a  study  done  on  this  tech¬ 
nique  to  establish  its  usefulness  in  comparison  with  an  exact  method.  From  this  work,  some 
insight  is  gained  into  the  extent  of  inhomogeneity  necessary  before  approximate  techniques  lose 
their  value.  It  is  now  widely  accepted  that  many  regions  of  the  Upper  Mantle  contain  presently 
active  or  fossil  slabs.  As  a  continuing  study  of  the  effect  these  inhomogeneities  have  on  seismic 
propagation,  we  report  on  the  effect  that  slabs  have  on  the  spectral  content  of  body-wave  signals. 
The  result  is  difficulty  in  distinguishing  slab  effects  from  anelastic  attenuation. 

Noise  studies  are  under  way  for  NORSAR  and  ALPA.  Long-period  noise  at  NORSAR  has 
the  expected  constituent  coming  from  the  Atlantic  Ocean,  but  also  a  significant  component  ap¬ 
pears  to  originate  to  the  east  of  NORSAR.  Such  noise,  which  is  in  the  same  sector  as  Eurasia, 
may  need  some  careful  attention  in  discriminating  small  events.  Joint  ALPA/LASA  noise  studies 
have  revealed  sources  of  ocean-generated  microseisms  which  are  close  to  regions  of  low  pres¬ 
sure  as  reported  on  weather  maps.  Long-period  signals  at  NORSAR  are  being  collected  for  the 
purpose  of  studying  variations  in  character  of  the  surface-wave  train  with  event  location.  In 
addition  to  these  variations,  it  has  been  possible  to  identify  portions  of  the  array  which  are  less 
sensitive  to  the  important  shorter-period  continental  phase.  First  array-processing  results 
are  given  for  NORSAR  discrimination  purposes.  Beamforming  and  matched-filter  techniques 
appear  to  be  working  well,  and  it  should  soon  be  possible  to  establish  a  surface-wave  threshold 
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Summary 


for  NORSAR;  this  figure  will  probably  be  seasonally  dependent.  LASA  corrections  for  core 
phase  arrivals  at  individual  subarrays  are  now  included  in  our  standard  Lincoln-based  programs, 
and  are  proving  valuable  in  improved  processing  capabilities. 

With  the  availability  of  NORSAR  and  ALPA  data,  it  is  desirable  to  consider  these  arrays 
not  only  in  isolation  but  in  terms  of  a  continent -wide  network  for  discrimination.  For  this  pur¬ 
pose,  two  new  projects  have  been  initiated.  A  search  is  being  carried  out  to  identify  the  most 
efficient  means  of  creating  a  bulletin  using  arrays  as  well  as  single  stations.  Thus  far,  use  of 
the  LASA  detection  log  in  conjunction  with  a  small  number  of  WWSSN  stations  is  proving  a  power¬ 
ful  means  of  obtaining  a  much  larger  earthquake  catalogue.  Capabilities  of  the  long-period 
constituent  of  the  WWSSN  for  discrimination  are  being  assessed  with  a  view  to  studying  the  larger 
population  we  expect  to  have  available  shortly. 

Recently,  we  have  been  receiving  data  from  the  Block-Moore  seismometer  installed  near 
San  Diego.  We  report  on  noise  and  signal  analysis  and  conclude  that  the  best  signal-to-noise 
ratio  at  this  site  is  in  the  0.03-  to  0.05 -Hz  band.  It  has  also  been  possible  to  generate  imitation 
WWSSN  records  from  these  data  in  order  to  obtain  an  understanding  of  the  degradation  that  too 
broad  a  filter  has  on  the  detectability  level  of  surface  waves.  The  relation  between  nonpropagating 
seismic  noise  and  microbarograph  output  has  been  studied  at  LASA.  We  concluded  that  correla¬ 
tion  techniques  enable  nonpropagating  noise  to  be  significantly  reduced. 

The  large  amount  of  data  now  available  in  our  Lincoln  facility  is  being  documented  for 
access  by  a  program  called  LISTAR.  This  program,  based  on  the  remote-console  accessed 
Lincoln  360,  should  enable  us  to  generate  lists  of  events  of  interest  within  minutes,  and  have 
full  information  on  how  to  find  relevant  data.  The  Analysis  Console  has  been  the  subject  of  a 
movie  designed  to  demonstrate  its  capabilities;  it  also  stars  in  a  BBC  television  production  on 
solid  earth  geophysics. 

Finally,  LASA  has  been  used  to  demonstrate  for  the  first  time  the  presence  of  seismic 
phase  PKJKP,  and  hence  the  rigidity  of  the  inner  core.  This  issue  has  been  in  doubt  for  35  years, 
but  it  now  seems  possible  to  assign  a  shear  velocity  of  2.95  ±0.1  km/sec  to  the  inner  core.  Fur¬ 
thermore,  LASA  is  being  used  to  study  fine  structure  on  the  boundary  of  the  inner  core. 

D.  Davies 
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GLOSSARY 


ALPA 

Alaskan  Long-Period  Array 

IGPP 

Institute  of  Geophysics  and  Planetary  Physics, 

University  of  California,  La  Jolla 

LASA 

Large  Aperture  Seismic  Array 

LP 

Long  Period 

LPE 

Long- Period  Experiment 

NORSAR 

Norwegian  Seismic  Array 

NOS 

National  Ocean  Survey 

SAAC 

Seismic  Array  Analysis  Center 

SATS 

Semiannual  Technical  Summary 

WWSSN 

World-Wide  Standard  Seismographic  Network 
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SEISMIC  DISCRIMINATION 


I.  THE  SEISMIC  SOURCE 

A.  OBSERVATIONS  OF  pP  IN  THE  SHORT-PERIOD  SOURCE  FUNCTION  OF  NTS 
EXPLOSIONS  RECORDED  AT  NORSAR 

At  teleseismie  distances,  the  effective  source  function  of  explosions  contains  not  only  the 

primary  P- phase,  but  also  the  pP-phase  whieh  is  reversed  in  polarity  and  arrives  within  a 

second  after  the  P-pulse  of  most  explosions.  A  third  pulse  is  often  seen  with  the  same  polarity 

as  the  P-pulse  and  later  than  pP.  The  detection  of  pP-phases  has  been  done  predominantly  by 

1-3 

frequency  domain  analyses  of  the  first  few  seconds  of  the  P-wave,  because  interference  be¬ 
tween  the  P  and  pP  makes  it  difficult  to  time  the  pP  onset.  Implicit  in  such  studies  is  the  assump¬ 
tion  that  the  pP-pulse  has  the  same  shape  as  the  P-pulse  exeept  for  the  change  of  polarity.  In 
the  frequency  domain,  this  produces  a  sharp  sealloping  of  the  amplitude  speetrum  of  the  wave¬ 
form,  giving  minima  at  periods  T^  =*  r/n  for  n  =  1,  2,  .  .  .  ,  where  t  is  the  delay  time  for  pP 

relative  to  P.  If,  however,  the  pP  has  a  different  shape  and  a  lower  amplitude  than  the  P-pulse 

3 

and  other  phases  are  present,  then  the  speetral  minima  may  disappear  or  be  inconclusive,  as 
is  often  the  ease. 

An  alternative  approach  to  the  separation  of  P-  from  pP-phases  is  to  simply  remove  the 
time-domain  instrument  response  from  the  observed  data  and  interpret  the  deconvolved  reeords 
as  displacement  source  functions.  In  this  manner,  one  obtains  a  direct  pieture  of  the  P,  pP  and 
other  phases  whieh  must  be  included  in  any  realistie  souree  function. 

This  is  being  attempted  first  with  explosions  of  known  depth  and  location.  Figure  1-1  shows 
a  suite  of  eleven  NTS  events  recorded  at  NORSAR.  The  shot  depths  indicated  for  these  events 

4 

are  taken  from  the  explosion  summary  by  Springer  and  Kinnaman.  The  waveforms  are  arranged 
in  order  of  increasing  depth  and  aligned  in  time  at  the  first  trough,  rather  than  at  the  first  peak. 
This  was  done  in  order  to  eliminate  the  effeet  of  the  variable  rise  time  of  the  first  motion,  and 
it  also  corresponds  to  the  maximum  displacement  of  the  impulse  displacement  response  of  the 
short- period  seismometer  at  NORSAR. 

Each  waveform  is  a  steered  beam  at  either  the  interim  array  of  NORSAR,  at  the  subarrays 
near  Faldalen  or  Oyer,  or  at  the  temporary  array  OONW  near  Lillehammer.  The  interim  array 
consists  of  a  single  seismometer  from  eaeh  subarray  of  the  complete  NORSAR  array. 

The  clearest  pP-phasc  oeeurs  on  the  seismogram  of  Stinger  at  Oyer.  Due  to  the  short  dura¬ 
tion  of  the  P-pulse,  the  pP  is  elearly  separated,  inverted  and  delayed  about  0.9  sec  relative  to 
P.  On  the  other  seismograms,  the  pP-phases  are  not  nearly  as  clear,  due  to  both  imperfeet 
reflection  of  P  at  the  free  surface  and  the  lower  frequency  content  of  the  primary  pulse.  On 
many  events  there  is  a  strong  seeond  trough,  which  has  the  wrong  polarity  and  time  for  pP. 

This  pulse  has  a  time  delay  which  increases  with  depth  and  is  part  of  the  effective  souree  function. 
For  example,  at  Oyer  there  is  increase  in  delay  time  of  almost  a  seeond  for  this  seeond  trough 
for  events  Purse  down  through  Roxear.  Similarly,  on  the  interim  array.  Hutch,  Rlenton  and 
Shaper  are  quite  similar,  whereas  Handley  has  a  seeond  trough  almost  a  seeond  later  than  the 
trough  in  the  three  shallower  shots. 
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Section  I 


The  seismograms  of  Fig.  1-1  were  deconvolved  by  removing  the  impulse  displacement 
response  of  the  short-period  seismograph  which  is  displayed  to  the  left  of  Stinger  and  Knox. 

As  shown,  this  response  has  been  convolved  with  an  attenuation  operator  with  t*  =  0.25  to  simu¬ 
late  most  of  the  average  earth  attenuation  from  the  NTS  area  to  NORSAR.5  The  deconvolved 
displacement  waveforms  are  shown  in  Fig.  1-2  aligned  in  time  at  the  first  peak.  Neglecting  the 
effect  of  layering  under  the  receiver,  these  seismograms  can  be  interpreted  as  the  effective 
far-field  displacement  source  for  each  event.  The  first  peak  of  each  wave  is  the  direct  P-pulse, 
followed  by  a  trough  which  is  the  pP-phase.  A  second  peak,  which  is  often  strong,  corresponds 
to  the  second  trough  seen  on  the  original  waveforms.  It  has  also  been  noted  by  Kulhanek  on 
other  NTS  shots  recorded  by  the  Swedish  seismic  network. 

Figure  1-3  shows  a  plot  of  the  delay  times  vs  depth  for  pP  picked  directly  from  the  decon¬ 
volved  seismograms.  The  receiver  array  is  indicated  for  each  event  to  show  that  the  pP  delay 
times  are  not  receiver-dependent.  Although  the  times  generally  increase  with  depth,  there  is 
considerable  scatter  in  the  data.  In  order  to  isolate  the  source-region  effects,  the  same  data  are 
plotted  in  Fig.  1-4,  with  symbols  indicating  the  NTS  source  region.  These  data  indicate  that 
Yucca  Flat  has  average  compressional  velocities  somewhat  slower  than  Pahute  Mesa. 

Average  compressional  velocities  for  the  media  above  each  shot  were  calculated  from  the 
formula 


n  =  2h/r 


where  h  is  the  shot  depth  in  kilometers,  and  t  is  the  observed  pP-P  delay.  The  data  for  these 
calculations  are  listed  in  Table  1-1.  Kulhanek  used  observed  spectral  minima  ~  1  Hz  to  infer 
pP  delay  times  for  NTS  shots  recorded  by  the  Swedish  seismic  network.  His  average  velocities 
and  the  velocities  m  Table  1-1  are  superposed  in  Fig.  1-5.  The  events  are  separated  into  three 
NTS  regions,  namely  Yucca  Flat,  Pahute  Mesa  and  Hot  Creek  Valley.  Of  the  two  sets  of  data, 
the  only  common  events  are  Boxcar  and  Faultless,  which  show  good  agreement  between  the  two 
methods. 


TABLE  1-1 

DELAY  TIME  t  VS  DEPTH  FOR  EVENTS 

STUDIED  AT  THREE  TEST  SITES 

Depth 

T 

* 

a 

NTS  Region 

Date 

Event 

(*) 

(sec) 

(km/sec) 

Pahute  Mesa 

5/9/69 

Purse 

599 

0.  75 

1.59 

3/22/68 

Stinger 

668 

0.85 

1.57 

4/26/68 

Boxcar 

1158 

1.  10 

2.  10 

3/26/70 

Ha  nd 1 ey 

1207 

1.  15 

2.09 

Yucca  Flat 

7/16/69 

Hutch 

549 

0.65 

1.69 

4/30/69 

Blentan 

557 

0.75 

1.49 

3/23/70 

Shaper 

561 

0.80 

1.40 

9/6/68 

Noggin 

582 

0.80 

1.46 

2/21/68 

Knox 

645 

0.95 

1.36 

10/18/67 

Lanpher 

715 

0.95 

1.51 

Hat  Creek 

Valley 

1/19/68 

Faultless 

975 

0.90 

2.  16 

*  Average  velocity  af  the  rocks  aver 

the  shat. 
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The  data  for  Pahute  Mesa  indicate  that  the  average  velocity  varies  linearly  from  about  1.5 
to  2. 5  km/sec  over  the  depth  range  600  to  1400  m.  Average  velocities  for  Yucca  Flat  vary  from 
about  1.1  to  1.7  km/sec  from  depths  of  460  to  760  m,  most  of  the  data  scattering  about  1.5  km/sec. 
These  velocities  are  reasonable  when  compared  with  a  cross  section  of  the  Yucca  Valley  given 
by  Springer/'  The  alluvium  surface  layer  varies  from  0.4  to  0.9  km/sec,  and  the  underlying  tuff 
stratum  has  velocities  up  to  1.5km/sec. 

Returning  to  Fig.  1-3,  it  appears  that  the  absence  of  a  clear  pP  on  the  original  seismograms 
can  be  explained  by  the  shape  of  the  pP-pulse  on  the  deconvolved  seismograms.  For  example, 
the  P-  and  pP-pulses  of  Stinger  appear  equally  sharp  in  this  figure  which  accounts  for  their 
clarity  on  the  original  seismogram.  However,  for  Boxcar  and  Handley,  the  pP-pulses  on  the 
deconvolved  seismograms  are  considerably  broader  and  of  lower  amplitude  than  the  P-pulses. 
This  imperfect  reflection  at  the  free  surface  produces  the  broad,  indistinct  pP-phases  seen  on 
the  original  seismograms. 

Finally,  we  note  that  the  seeond  peak  on  the  deconvolved  seismograms  is  often  large,  for 
example  on  Rlcnton  and  Shaper.  This  pulse  may  be  due  to  the  spalling  of  the  ground  surface  over 
the  shot  and  should  be  included  as  part  of  the  effective  source  function,  observed  at  teleseismie 
distances.  An  example  of  spalling  is  shown  in  Fig.  1-6.  This  shows  the  vertical  components  of 
acceleration  and  displacement  observed  on  the  ground  surface  over  a  "typical”  test  shot/  which 
unfortunately  is  not  identified.  The  first  aeeeleration  peak  is  the  direct  wave  arriving  at  the 
surface,  which  after  reflection  produces  the  pP-phase  observed  telcseismically.  The  second 
peak  is  caused  by  the  "slap  down"  of  the  rock  material  after  about  a  second  of  free  fall.  This 
acts  like  a  secondary  source  which  then  radiates  out  telcseismieally  with  the  same  polarity  as 
the  primary  P-pulsc.  c.  w>  Frasier 


B.  SOURCE  FUNCTION  FOR  EXPLOSIONS  FROM  CORE  REFLECTIONS 


In  the  last  two  SATS  (31  December  1970,  DDC  AD-718971;  and  30  June  1971,  L)DC  AD- 
728210),  I  discussed  in  some  detail  the  problem  of  looking  at  source-time  functions  for  explo¬ 
sions.  This  problem  is  important  for  understanding  the  physical  basis  for  discrimination.  It 
has  been  elear  for  some  time  that  the  teleseismie  P-wave,  although  a  clear  and  simple  signal, 
often  is  overprinted  with  confusing  detail  of  lateral  heterogeneity  in  the  mantle.  The  best  way 
to  look  seismically  at  the  source  seems  to  be  by  means  of  a  core  reflection,  for  which  P-waves 
have  emerged  nearly  vertically  and  are  thus  least  prone  to  see  lateral  variations.  In  earlier 
reports,  the  phase  PKiKP,  reflected  from  the  inner  core,  has  been  used.  The  assumption  was 
made  that  the  inner-core  boundary  was  sharp  and,  thus,  that  reflection  from  it  gave  as  undis¬ 
torted  a  picture  of  the  short-period  source  function  as  seismology  could  provide.  Recent  work 
on  Milrow  and  Cannikin  has  shown  this  assumption  to  be  incorrect.  As  a  result  of  new  data, 
one  ean  now  construct  both  a  revised  source  function  and  a  tentative  structure  for  the  inner-eore 
boundary  (see  Sec.  VII). 

Figure  I-7(a-b)  shows  LASA  signals  for  phases  PeP,  PKiKP  and  PKKKP  from  Milrow  and 
Cannikin.  PKKKP  is  particularly  remarkable.  It  has  a  dT/dA  of  2.3  sec/dcg  and  a  travel  time 

of  41m  51s  and  appears  to  belong  near  the  II  point  of  core  travel-time  tables  (A  =  408°).  As 

2  2 

far  as  I  am  aware,  this  is  the  first  report  of  PKKKP^j.  In  this  region  d  T/dA  is  very  high, 
which  explains  why  the  phase  is  visible  at  all. 
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The  first  thing  to  be  said  about  these  data  is  that  PKKKP  and  PcP  are  both  very  short  signals 

lasting  less  than  2  sec.  PKKKP  has  a  reversed  polarity  which  possibly  arises  from  double  pas- 

q 

sage  through  a  caustic  in  the  core  (a  minimax  reflection  as  described  by  Jeffreys  and  Lapwood  ). 
This  point  needs  further  investigation.  Nevertheless,  it  seems  reasonable  to  assert  that  PKKKP 
and  PcP  represent  a  better  view  of  the  source  than  PKiKP.  We  can  now  compare  these  two 
phases  for  Milrow  and  Cannikin.  No  yields  have  been  published  for  these  events, but  they  are 
presumed  to  have  yields  of  1  and  5  megatons,  respectively.  It  is  clear  from  Fig.  I-7(a-b)  that 
the  time  scales  as  well  as  the  amplitude  scales  of  the  two  PKKKP  signals  are  different.  How¬ 
ever,  the  scaling  is  very  much  less  than  5  to  1.  The  depths  of  burial  of  the  two  events  were 
4000  feet  (Milrow)  and  6000  feet  (Cannikin)  and  the  time  scaling  is  roughly  in  this  ratio,  sug¬ 
gesting  that  depth  of  burial  may  have  an  effect.  Unfortunately,  it  is  difficult  to  sort  out  scaling 
owing  to  burial  depth  from  scaling  simply  arising  from  explosion  size,  which  will  go  as  the  cube 
root  of  yield.  However,  it  is  fairly  clear  from  Fig.  1-7  that  the  PcP  and  PKKKP  signals  both 
comprise  more  than  one  discrete  pulse.  Deconvolving  the  impulse  displacement  response  with 
attenuation  from  the  PcP  signal  yields  a  tentative  source  model  for  Milrow,  shown  in  Fig.  l-8(a-b). 

Compare  Fig.  1-8  with  the  results  of  Frasier  in  Sec.  A  above.  The  agreement  is  remark¬ 
ably  close.  We  each  conclude  that  the  signal  radiated  to  teleseismic  distances  consists  of  a 
primary  pulse  together  with  a  complicated  combination  of  a  reflected  pulse  (negative)  and  a  signal 
contributed  by  spalling  (positive).  The  ground  displacement  beneath  the  shot  thus  appears  from 
teleseismic  data  to  contain  an  impulsive  term  together  with  a  step  function  (not  necessarily  of 
the  same  size). 

I  wish  to  acknowledge  Dr.  E.  A.  Flinn  (Teledyne-Geotech)  for  providing  the  Cannikin  data. 

D.  Davies 


C.  SIZE  OF  EARTHQUAKES 


The  dimensions  of  the  rupture  region  of  an  earthquake  are  difficult  to  acquire  from  seismic 
data.  Even  if  the  earthquake  breaks  surface  in  an  accessible  region,  this  gives  no  clear  indica¬ 
tion  of  whether  the  visible  rupture  length  is  the  total  length  of  activity  and,  of  course,  no  guide 
at  all  to  the  depth  of  the  ruptured  region.  Study  of  the  area  of  aftershock  regions  has  been  used 
for  large  earthquakes.10  Indirect  methods  have  become  popular  in  the  last  year  or  two,  depend¬ 
ing  mainly  on  spectral  observations  of  seismic  signals,  and  the  determination  of  a  11  corner 
frequency”  which  is  related  to  the  size  of  an  earthquake.  We  propose  a  simple  time-domain 
approach  to  the  problem  which  so  far  has  yielded  interesting  results. 

It  is  assumed  that  an  earthquake  can  be  described  as  a  plane  region  across  which  sudden 
rupture  occurs.  This  rupture  nucleates  at  a  certain  point  on  the  fault  plane  and  spreads  over 
the  plane  until  it  is  terminated  by  some  unspecified  obstacle.  The  rupture  is  often  presumed 
to  propagate  with  a  velocity  less  than  that  of  the  shear  wave  velocity  of  the  medium.  Little  is 
known  about  the  "quenching”  process  by  which  the  rupture  is  brought  to  a  standstill.  If  this 
description  of  an  earthquake  is  valid,  the  duration  (t)  of  the  far-field  radiated  signal  should 
bear  a  close  relation  to  the  size  of  the  event.  In  particular,  one  would  expect  that  by  observa¬ 
tion  with  P-waves  having  near-source  velocity  a 
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where  L  is  the  largest  distance  the  rupture  propagates  with  velocity  v^.  The  extreme  (equal) 
cases  correspond  to  a  rupture  propagating  in  the  same  or  opposite  direction  as  the  P-wave.  A 
reasonable  upper  bound  to  the  largest  dimension  of  the  earthquake  is  then  2L,  or  typically  16  t 
kilometers.  A  more  probable  value  if  fault  plane  and  rupture  velocity  vector  orientation  are 
random  is  8  r  kilometers.  In  the  past,  determination  of  r  (which  clearly  needs  short-period 
seismograms)  has  been  severely  limited  by  the  difficulties  associated  with  short-period  data. 
However,  a  digital  array  removes  many  of  these  difficulties  and  high-quality  records  are  readily 
available.  One  problem  remains  —  that  lateral  heterogeneities  in  the  upper  mantle,  and  partic¬ 
ularly  in  the  source  region,  lead  to  very  confusing  coda  and  do  not  allow  a  judgment  to  be  made 
of  where  the  signal  ends  and  the  coda  begins.  Here,  the  power  of  an  array  to  beamform  toward 
specified  phases  of  an  event  can  be  used  to  advantage.  Once  an  event  has  been  chosen,  several 
phases  can  be  examined  and  the  shortest  in  duration  selected  to  give  an  estimate  of  r.  Wc  find 
that  core  reflections  are  invariably  the  most  useful  for  this  purpose  as  the  initial  signal  has 
emerged  steeply  from  the  heterogeneous  zone. 

In  Fig.  I-9(a-o),  we  give  three  examples;  in  each  ease,  we  show  the  P-wave  and  PcP-wave 
beams  formed  by  LASA.  In  some  instances,  a  surface  reflection  is  also  visible.  It  is  clear  that 
in  no  case  is  the  signal  in  the  ground  at  LASA  of  more  than  1  sec  in  duration,  from  which  we 
infer  that  the  sourees  are  at  most  16  km  in  size  and  probably  much  smaller.  Particularly 
striking  is  event  (c).  The  signal  looks  exactly  like  the  instrument  response  to  a  step  function  in 
displacement,  and  we  propose  that  the  dimensions  of  this  deep  earthquake  are  at  most  4  km. 

The  difference  in  character  of  PcP  among  the  three  events  is  remarkable.  It  bears  infor¬ 
mation  on  the  time-domain  response  of  the  core-mantle  boundary,  and  we  are  constructing  a 

model  consistent  with  these  data.  .. 

I).  Davies 

A.  M.  Ziolkowski 
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Fig.  1-1  .  Eleven  NTS  P-waves  recorded  at  arrays 
in  Norway.  Events  are  arranged  by  increasing  depth 
in  meters  from  top  to  bottom.  To  left  afdata  is  shown 
shart-periad  displacement  response  af  seismometer, 
attenuated  ta  approximate  Earth's  absorption  fram 
NTS  to  Norwoy. 


Fig.  1-2.  Deconvolved  displacement  recards 
af  seismograms  in  Fig.  1-1 .  These  waveforms 
are  obtained  by  removing  instrument  response 
fram  data.  First  troughs  are  interpreted  ta  be 
pP-phases. 
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Fig.  1-3.  pP  delay  times  vs  depth 
for  NTS  events.  These  are  obtained 
by  inspection  af  Fig.  1-2.  Array  re¬ 
cording  each  event  is  shown. 


Fig.  1-4.  pP  delay  times  vs  depth 
far  NTS  events.  NTS  source  region 
is  indicated  far  each  event. 
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Fig.  1-5.  Comparison  af  average  velocities  at  NTS  regions  obtained  fram  pP  delay 
times.  Dats  obtained  by  Kulhanek  fram  observed  minima  in  amplitude  spectra  af 
NTS  shots.  Open  circles  calculated  using  delay  times  from  reduced  records  af  Fig.  1-2. 
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Fig.  1-6.  Surface  acceleration  and  displacement 
aver  a  “typical"  contained  NTS  explosion.  Sec¬ 
ond  acceleration  peak  occurs  when  rack  material 
returns  ta  ground  level  after  sec  af  free  fall. 
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Fig.  1-7.  Core  phases  from  (a)  Milraw  and  (b)  Cannikin. 
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Fig.  1-9.  P-  and  PcP-wave  beams  farmed  by  LASA  fram  three  earthquakes,  (a)  Kamondarsky  Islands, 
20  Januory  1  969;  h  =  23  km ,  m^  =  6. 1 .  (b)  Komandarsky  Islands,  4  April  1  969;  h  =  27  km,  mb  =  5.4. 
(c)  Alaska,  1 0  June  1 968;  h  =  1  82  km,  =  5.6. 
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II.  PROPAGATION 


A.  NUMERICAL  STUDIES  OF  SEISMIC  WAVE  PROPAGATION 

A  numerical  integration  technique  has  been  developed^  which  solves  the  elastic  equations 
governing  coupled  compressional  and  shear  waves  in  inhomogeneous  media.  By  neglecting  doubly 
backscattered  waves,  the  equations  reduce  to  a  set  of  first-order  differential  equations  which 

are  easily  numerically  integrated  for  media  in  which  standard  numerical  techniques  cannot  be 

2 

applied.  For  waves  passing  through  a  medium,  the  theory  is  valid  provided  that  {oo/a')  »  1 
2 

and  (u>//B)  »  1,  where  u>  is  the  angular  frequency  and  cv*  and  /3 1  are  the  gradients  of  the  com¬ 
pressional  and  shear  velocities,  respectively.  The  limits  to  which  these  high-frequency  re¬ 
strictions  can  be  pushed  are  examined  below. 

In  studies  of  seismic  source  functions,  it  is  important  to  consider  the  effect  of  velocity  and 
density  transition  zones  in  the  Earth  on  short-period  signals.  As  a  test  of  the  above  technique, 
the  plane- wave  transmission  response  for  a  transition  zone  was  calculated  and  compared  with 
the  discrete  time  solution,^  which  has  been  shown  to  converge  to  the  exact  continuous  time  so¬ 
lution.  The  transition  layer  is  shown  in  Fig.  II- 1.  Both  velocity  &  and  density  p  decrease 
linearly  over  a  2- km  zone,  and  have  end  values  which  are  typical  for  a  crustal  sedimentary  layer 
over  basement  rock.  Figure  II—  2  shows  the  discrete  time  transmission  response  of  the  transition 
to  an  incident  delta  function  of  compressional  particle  velocity  with  energy  density  1.  Three 
solutions  are  shown  for  discrete  time  intervals  At  =  0.4,  0.2,  and  0.1,  and  are  nearly  identical 
at  coincident  sample  points.  The  first  arrival  is  a  delta  function  of  area  ~0.99,  followed  by  a 
smooth  finite  amplitude  tail  which  is  the  accumulation  of  doubly  backscattered  multiples  which 
travel  forward  and  arrive  after  the  first  pulse.  Such  backscattering  is  neglected  in  the  numerical 
integration  scheme. 

Figure  11-3  shows  the  amplitude  and  phase  response  of  the  transmitted  wave  calculated  by 
each  method  as  a  function  of  frequency.  The  crosses  indicate  the  superposed  discrete  time 
solutions  (all  three  are  identical),  and  the  circles  show  the  numerical  integration  solutions. 

From  the  discrete  solution,  we  see  that  the  transmitted  delta  function  makes  the  spectrum  nearly 
flat  except  below  1  IIz,  where  the  negative  portion  of  the  response  in  Fig.  II- 2  lowers  the  ampli¬ 
tude.  The  phase  response  is  essentially  0.  This  verifies  the  assumption  in  the  numerical  inte¬ 
gration  scheme  that  backscattered  reflections  are  practically  negligible.  For  the  transition 
model  considered,  aj/cv’  =  10  corresponds  to  a  frequency  of  ~1.6llz.  Above  this  frequency, 

(cc/a')  >  1  00  »  1  and  the  numerical  integration  scheme  gives  results  which  closely  match  the 

discrete  time  solution.  Below  1.6  IIz,  the  amplitude  response  of  the  numerical  integration  scheme 
falls  off  rapidly,  as  it  should;  however,  the  phase  distortion  remains  less  than  30°.  This  implies 
that  time-domain  seismograms  can  be  synthesized  from  Fourier  components  of  the  transmission 
response  without  much  distortion  of  the  wave  shape. 


t  T.  E.  Landers,  "Elastic  Waves  in  Laterally  Inhomogeneous  Media,"  Ph.  D.  Dissertation, 
Stanford  University  (1971). 

}C.  W.  Frasier,  "Discrete  Time  Solution  of  Plane  P-SV  Waves  in  a  Plane  Layered  Medium," 
Geophysics  35,  197-219  (1970). 
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For  geologic  zones  with  less  severe  velocity  gradients,  such  as  the  Moho,  the  numerical 
integration  will  be  accurate  across  the  complete  short-period  band. 

Similar  studies  are  being  made  on  the  coupled  conversions  of  P-  to  S- waves  in  transition 
zones  for  both  the  transmission  and  reflection  problems. 

C.  W.  Frasier 
T.  E.  Landers 


B.  SPECTRAL  DISTORTION  WITHIN  SHADOWS 

The  effects  upon  seismic  wave  propagation  of  lateral  variations  in  the  Earth’s  structure 
have  been  investigated  in  several  previous  SATS  using  first-order  ray  theory,  and  it  was  found 
that  for  realistic  models  quite  large  amplitude  anomalies  and  even  shadows  can  be  produced. 
However,  ray  theory  is  strictly  applicable  only  in  the  limit,  as  the  frequency  tends  to  infinity. 

This  means  that  the  amplitude  anomalies  will  be  functions  of  frequency  or,  stated  differently, 
that  lateral  variations  can  distort  the  spectra  of  seismic  signals.  Here  we  present  some  esti¬ 
mates  of  the  magnitude  of  such  effects  for  the  case  of  a  dipping  high-velocity  lithospheric  slab 
beneath  an  island  arc. 

The  method  used  consists  of  orienting  the  coordinate  axes  so  that  the  velocity  in  the  vicinity 
of  the  slab  may  be  approximated  as  a  function  of  one  coordinate  only  (see  Fig.  II- 4).  The  problem 
is  thus  reduced  to  one  of  reflection  and  transmission  from  a  one- dimensional  velocity  variation, 
which  is  solved  by  direct  numerical  integration  of  the  equations  of  motion  resulting  when  the 
factor  exp[i(o)t  -  kx)]  is  removed. 

Here  we  consider  a  symmetric  model  in  which  the  two  seismic  velocities  and  the  density 

2 

are  given  by  functions  of  the  form  A  +  B  exp[-(z/h)  ].  Such  a  restriction  is  not  necessary;  any 
variation  at  all  can  be  handled  by  the  method.  Figure  11-5  shows  the  calculated  transmission 
coefficient  as  a  function  of  wavelength  for  P-waves  incident  at  a  variety  of  angles  upon  a  slab 
with  a  10- percent  velocity  increase.  Hay  theory  predicts  perfect  transmission  for  angles  less 
than  64.4°,  and  total  reflection  for  greater  angles;  thus,  all  the  cases  shown  lie  within  the 
geometrical  shadow.  It  is  clear  that  the  frequency  above  which  ray  theory  applies  is  very  roughly 
given  by  \  =  h,  but  is  a  strongly  varying  function  of  the  angle  of  incidence,  particularly  near 
grazing  incidence.  It  is  also  clear  that  in  a  shadow  zone  a  pronounced  tilt  will  be  introduced  into 
the  spectrum, which  will  be  difficult  to  distinguish  from  the  effect  of  anelastic  attenuation.  (The 
effect  of  attenuation  on  a  plot  such  as  Fig.  11-5  would  be  a  straight  line  passing  through  the 
upper  left-hand  corner.)  The  curve  for  i  =  71.0°  in  Fig.  11-5  corresponds  roughly  to  a  value  of  t* 
(T/Q,  where  T  is  travel  time  and  Q  is  the  effective  quality  factor)  of  about  6,  for  example. 

Thus,  not  only  magnitude  measurements  but  also  spectral  measurements  will  be  badly  biased 
when  shadowing  occurs.  This  effect  may,  in  fact,  prove  of  value  for  determining  the  thickness 
of  the  descending  lithosphere,  current  estimates  run  all  the  way  from  20  to  more  than  100  km. 

B.  R.  Julian 
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Fig.ll-1.  Tronsition  (oyer  between  two  ho  If  spaces.  Incident  P-wave  iso  delta  function 
in  higher-velocity  half  space.  Transmitted  wave  emerges  into  lower-velocity  half  space. 


Fig.  11-2.  Discrete  time  solution  for  transmitted  wove.  Three  different  time  increments  At  give 
neorly  identicol  volues  at  coincident  times.  First  spike  hos  oreo  of  ~0.99,  andnegotive  toil  hos 
areo  of  ~0. 06. 
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Fig.  11-3.  Comparison  af  transmission  responses  calculated  by  two  methods.  Circles  are  values 
obtained  by  numerical  integration,  and  crosses  shaw  discrete  time  solution.  Far  frequenc  iesabave 
u/a1  =  10,  numerical  integration  method  is  quite  accurate. 
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Fig.  11-4.  Schematic  diagram  af  descending  lithospheric  slab,  showing  orientation 
af  coordinate  system. 


ft/  X 

Fig.  11-5.  P-wave  transmission  coefficient  far  slab  model  with  10-percent  increase  in  campressianal 
and  shear  velocity  and  na  density  change,  as  a  function  af  wavelength  af  incident  wave  far  various 
angles  af  incidence.  Critical  angle  in  this  case  is  64.4°;  h  is  the  half-width  af  the  slab  (see  text). 
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A.  ANALYSIS  OF  LONG-PERIOD  MICROSEISMIC  NOISE  AT  NORSAR 

The  analysis  of  microseismic  noise  is  extremely  important  for  the  problem  of  seismic 
discrimination.  This  analysis  provides  the  data  to  determine  the  detection  capability  of  a  seis¬ 
mic  station,  or  array,  for  both  short-  and  long-period  signals.  In  addition,  the  analysis  provides 
information  concerning  the  character  of  the  seismie  noise.  This  information,  in  turn,  can  be 
used  to  determine  the  effectiveness  of  various  signal-processing  methods  in  reducing  the  noise 
power. 

There  have  been  several  investigations^  ^  of  the  seismie  noise  field  at  LASA.  These  anal¬ 
yses  have  lent  considerable  insight  into  the  nature  of  the  microseismic  noise  field.  We  will 

now  present  the  results  of  analysis  of  the  microseismic  noise  in  the  6-  to  40-sec  period  range 

4 

at  NORSAR.  High-resolution  frequency-wavenumber  spectra  have  been  measured  for  the  noise, 
as  well  as  power  spectra  and  coherence. 

The  frcqucncy-wavenumber  spectra  of  the  long-period  microseismic  noise,  in  the  6-  to 
40-sec  period  range,  have  been  measured  at  the  long-period  array  at  NORSAR.  This  array  con¬ 
sists  of  22  sets  of  3-component  long-period  seismometers  loeated  within  an  aperture  of  about 
125  km.  The  beam  pattern  of  the  array  in  k-space  is  shown  in  Fig.  Ill— 1 .  The  high-resolution 

method  was  employed,  and  the  parameters  of  the  measurement  were  the  same  as  those  used 

4 

previously  for  the  analysis  of  long-period  noise.  Five  different  noise  samples  were  analyzed 
during  a  1 -month  interval  in  June  1971.  The  results  to  be  presented  are  typical  of  the  entire 
set  of  data. 

The  results  of  one  of  the  measurements  for  a  1  June  1971  noise  sample  are  shown  in 
Fig.  IIl-2(a-c)  for  frequencies  of  0.03,  0.04  and  0.05Hz,  respectively,  as  measured  for  the  ver¬ 
tical  instruments.  The  corresponding  data  for  the  north-south  and  east-west  seismometers  are 
shown  in  Figs.  III-3(a-c)  and  -4{a-c),  respectively.  A  typical  power  spectrum  is  shown  in 
Fig.  1II-5.  The  method  used  to  estimate  the  eoherenee  and  spectra  is  the  same  as  that  given 
previously.^  The  data  in  Fig.  III-2(a-c)  indicate  that  the  20-  to  40-sec  microseisms  contain  a 
component  which  consists  of  fundamental-mode  Rayleigh  waves  whieh  propagate  from  a  wide 
range  of  azimuths.  This  fact  is  also  corroborated  by  the  data  in  Fig.  III-3(a-e).  The  data  in 
Fig.  III-4(c)  show  that  the  noise  also  consists  of  fundamental-mode  Love  waves  which  arrive  at 
the  array  from  the  same  direction  as  the  fundamental-mode  Rayleigh  waves.  This  indicates 
that  the  same  mechanism  may  be  responsible  for  the  generation  of  both  Rayleigh-  and  lave-wave 
microseismic  noise  at  NORSAR. 

The  coherence  data  show'  that  the  20-  to  40-sec  microseisms  also  contain  a  nonpropagating 
component.  The  noise  levels  depicted  by  the  power  speetra  in  Fig.  Ill—  5  are  typical  of  those 
found  previously  at  IASA.^  j  £ap0n 


B.  SIMULTANEOUS  ANALYSIS  OF  LONG-PERIOD  MICROSEISMIC  NOISE 
AT  ALPA  AND  LASA 

The  long-period  background  seismic  noise  limits  the  identification  threshold  at  which  the 
powerful  Mg-m^  discriminant  can  be  applied  at  LASA.  This  noise  has  been  investigated^  and 
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found  to  eonsist  of  two  components.  One  component  propagates  across  LASA  as  a  fundamental¬ 
mode  Rayleigh  wave  and  is  believed  to  be  caused  by  the  action  of  surf  on  coastlines.  The  other 
component  is  nonpropagating  (or  incoherent)  noise,  and  data  have  been  obtained  that  indicate  it 
is  eaused  by  atmospheric  pressure  fluctuations.  This  component  always  provides  a  significant 
contribution  to  the  total  noise  power. 

The  previous  analyses  were  able  to  obtain  only  the  direction  of  arrival  at  IJ\SA  of  the  prop¬ 
agating  component  of  the  noise.  Thus,  it  was  not  possible  to  determine  the  origin  of  these  sur¬ 
face  waves.  Such  a  determination  would  be  possible  by  obtaining  the  direction  of  arrival  of  these 
surface  waves  simultaneously  at  two  arrays  located  reasonably  elose  to  eaeh  other.  This  meas¬ 
urement  has  become  feasible  to  perform  as  a  eonsequence  of  the  availability  of  data  from  the 
Alaskan  Long -Period  Array  (ALPA). 

The  ALPA  is  an  array  of  nineteen  3-eomponent  long-period  seismometers  arranged  20  km 
apart  in  a  filled  hexagonal  pattern  of  about  80-km  diameter.  The  eenter  of  the  array  is  loeated 
about  50km  north  of  Fairbanks,  Alaska.  The  beam  pattern  of  the  ALPA  is  depieted  in  Fig.  1II-6 
where  it  is  seen  that  there  is  a  serious  aliasing  sidelobe  whose  magnitude  is  only  1  dB  less  than 
that  of  the  mainlobe.  This  severe  sidelobe  problem  is  due  to  the  almost  perfect  periodieity  in 
spaee  of  the  seismometers  arranged  in  the  filled  hexagonal  pattern.  This  sidelobe  is  loeated  at 
a  distance,  in  wavenumber  spaee,  from  the  mainlobe  of  about  0.057  Ilz/km.  Thus,  there  eould 
be  serious  problems  in  beamforming  or  frequeney-wavenumber  analysis  of  surfaee  waves  with 
a  period  of  l/0.5  X  0.057  X  3.7)  =  10  see.  This  eould  be  a  serious  problem  for  seismie  discrim¬ 
ination  using  the  M^-m^  criterion,  sinee  Rayleigh  waves  observed  at  short  distances  from  their 
origin  often  have  their  dominant  energy  at  the  shorter  periods,  sueh  as  in  the  10-  to  14-sec 
period  range. 

4 

The  high-resolution  frequency-wavenumber  speetra,  eoherenee  and  power  speetra  of  the 

noise  in  the  6-  to  40-sec  period  range  were  measured  simultaneously  at  ALPA  and  LASA  for 

time  periods  during  the  days  16  September,  25  September  and  20  Oetober  1970.  The  results  for 

two  of  these  time  periods  are  shown  in  Figs.  Ill— 7  and  -8.  The  directions  shown  in  these  figures 

pertain  to  seismie  noise  in  the  20-  to  2  5-sec  period  range  and  were  obtained  from  the  frequency- 

wavenumber  spectra  results  for  the  vertical,  north-south  and  east -west  arrays  at  both  ALPA 

and  LASA.  In  Figs.  Ill— 7  and  -8,  we  see  that  the  location  of  the  noise  source  extends  over  a  wide 

region  of  the  earth.  In  addition,  this  region  contains  a  low-atmospherie-pressurc  area.  The 

locations  of  these  areas  were  obtained  from  surface  eharts  showing  isobars  and  low-  and  high- 

pressure  regions.  Thus,  the  results  of  Figs.  Ill— 7  and  -8  show  that  the  souree  of  the  seismic 

noise  ean  be  associated  with  the  existence  of  low-pressure  atmospheric  disturbances  near  eoast- 

2 

lines,  as  indicated  by  Haubrieh  and  McCamy.  The  results  for  the  25  September  1970  noise  sam¬ 
ple  did  not  yield  a  location  of  the  noise  souree  whieh  eould  be  associated  with  a  low-pressure 
region. 

The  power  speetra  of  the  16  September  197  0  noise,  for  the  vertical,  north-south  and  cast- 
west  seismometers  at  site  2-4  at  ALPA  and  site  A0  at  LASA,  are  shown  in  Figs.  III-9(a)  and  (b), 
respectively.  These  data  show  that  the  seismie  noise  levels  at  ALPA  and  LASA  are  comparable. 
Similar  results  were  obtained  for  the  other  two  noise  samples  that  were  analyzed. 

J.  Capon 
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C.  DISCUSSION  OF  A  SUITE  OF  LONG-PERIOD  SIGNALS  RECORDED 
IN  NORWAY  FROM  ASIAN  EVENTS 

Lincoln  has  reeeived  a  limited  set  of  long-period  digital  recordings  from  the  completed 
NORSAR  and  a  set  of  plots  of  vertical  instruments  for  four  events  from  Central  Asia.  Only  one 
well-reeorded  event  from  the  area  has  been  found  on  the  digital  tapes.  Data  from  these  five 
events  have  been  combined  with  recordings  of  thirteen  additional  1968  events  obtained  during  the 
noise-survey  stage  of  NORSAR  to  make  some  preliminary  observations  eoneerning  the  events, 
the  site  and  Mg-m^  discrimination  procedures.  Figure  111-10  shows  the  location  of  the  eighteen 
events  in  this  study  and  Table  III— 1  gives  more  detailed  information  eoneerning  the  individual 
events. 

The  1968  data  eame  from  sites  identified  as  Oyer,  Trysl  and  Faldalen.  The  Oyer  site  has 
now  beeome  subarray  01 C.  The  Faldalen  site  was  near  subarray  07C,  and  Trysl  was  outside 
the  eurrent  NORSAR  to  the  east  and  slightly  north.  The  current  NORSAR  geometry  and  site 
nomenclature  are  shown  in  Fig.  Ill  —  1 1 . 

Figure  III— 1  2  shows  long-period  vertical  Rayleigh  waves  recorded  at  all  sites  in  NORSAR. 
Earthquake  8  and  presumed  explosion  3  show  a  relatively  well-developed  phase  with  a  period  on 
the  reeord  of  about  14  see.  This  signal  is  the  burst  of  higher  frequency  which  is  followed  by 
lower-frequeney  arrivals.  We  note  in  passing  the  strong  similarity  of  signal  shape  between  these 
two  events.  Since  the  paths  are  similar,  this  suggests  long-period  excitation  with  similar-shaped 
speetra  for  these  two  events.  Event  10  has  only  the  slightest  hint  of  this  14-see  signal,  although 
the  path  is  again  not  very  different.  Though  this  eould  be  due  to  path  effects,  it  is  more  likely 
that  the  source  is  rieher  in  the  lower  frequencies  so  the  14-see  signal  has  a  relatively  smaller 
contribution.  There  is  no  hint  of  the  14-see  signal  for  event  2.  The  burst  of  largest  signal  is 
at  about  an  18-see  period  and  may  correspond  to  a  flat  portion  of  the  group  velocity  curve  (see 
Filson’s  comments  on  group  velocity  in  See.  D  below). 

Event  3  shows  a  very  interesting  fact  about  the  14-see  phase.  Specifically,  note  that  its 
contribution,  relative  to  other  periods,  is  much  less  at  sites  01 C  and  02C  than  at  many  others. 

In  fact,  there  is  a  strong  tendency  for  it  to  be  weak  in  the  northwest  part  of  the  array.  The  ef¬ 
fect  is  strong  for  a  relatively  small  change  in  rceeiver  site. 

Figure  II I  —  1  3  shows  Rayleigh  and  Love  signals  for  the  remaining  events.  The  vertical, 
north-south,  and  east-west  components  are  shown  for  eaeh  event.  Note  that  data  were  obtained 
only  from  the  Oyer  site  for  many  of  these  events  so  that  the  high-frequeney  component  is  expected 
to  be  relatively  small  or  shifted  slightly  to  longer  periods. 

The  data  have  been  analyzed  to  obtain  measurements  of  surfaee-wave  magnitude  using  both 
Rayleigh  and  Love  waves.  The  results  are  tabulated  in  Table  III— 1 .  Figure  lll-14(a)  is  a  plot 
of  Rayleigh-wave  Mg  vs  the  NOS  m^  values.  Mg  was  computed  from  Mg  =  log(A/T)  +  1.66  log  A, 
where  A  is  peak-to-peak  signal  in  millimicrons,  A  is  distance  in  degrees,  and  T  is  period  in 
seconds.  The  period  was  restricted  to  the  18-  to  23-see  band.  The  same  formula  was  used  for 
Love  waves  but  the  largest  signal  with  T  >  18  was  used.  These  data  are  shown  in  Fig.  111—14 (b) . 

Ix)ve  waves  reeeived  from  explosions  at  a  receiver  are  usually  smaller  than  the  Rayleigh 
waves,  but  there  is  a  reasonable  number  of  earthquakes  for  whieh  this  is  not  true.  Thus,  if  the 
larger  of  the  Rayleigh  and  L,ove  Mg  values  is  plotted  vs  m^,  the  separation  should  be  improved 
over  a  customary  IU  -m^  plot  such  as  Fig.  lll-14(a).  This  has  been  done  in  Fig.  III-14(c)  and, 
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although  the  population  is  limited,  we  do  see  some  improvement.  The  procedure  makes  sense 
theoretically  as  well  as  experimentally  since  the  Rayleigh  and  Love  radiation  patterns  from 
earthquakes  are  45°  out  of  phase,  thus  a  strong  Love  wave  may  be  received  with  a  weak  Rayleigh 
wave  even  if  the  aetual  generation  of  Rayleigh  and  Love  waves  is  not  anomalous.  Since  horizon¬ 
tal  noise  at  NORSAR  is  often  as  small  as  the  vertical,  this  use  of  Love  waves  is  quite  promising. 

Two  events  —  one  presumed  explosion  and  one  earthquake  —  could  be  problems  with  these 
Mg-m^  data.  These  are  events  13  and  15  and  neither  was  affected  by  taking  the  larger  of  the 

Love  and  Rayleigh  M  values.  Both  of  these  events,  as  well  as  five  other  events  of  those  with 

s  5 

M^-m^  values  on  Fig.  III-14,  were  included  in  a  previous  IASA  experiment.  values  or  bounds 

from  that  experiment  are  also  included  in  Table  III—  1 .  Figure  III -1 4(d)  is  an  Mg-m^  plot  using 
the  largest  of  the  two  NORSAR  values  or  the  IASA  value  for  M^.  The  plot  is  unchanged  from 
Fig.  Ill— 1 4(c)  except  that  earthquake  15  has  been  raised  slightly.  The  improvement  is  not  large 
but  encourages  us  to  pursue  such  extremum  methods  more  thoroughly.  Our  motivation  in  taking 
the  maximum  over  stations  is  that  earthquakes  should  have  a  more  azimuthally  variable  radiation 
pattern  than  explosions,  and  separation  may  thus  be  improved  over  using  simple  averages.  Of 
course,  using  maxima  will  require  careful  application  of  path  corrections  and  learning  system¬ 
atic  effects  other  than  radiation  patterns  which  could  cause  a  station  to  have  a  high  value  of  Mg. 
This  has  not  really  been  done  here  in  combining  the  IASA  and  Norway  data.  The  point  is  to  de-. 
velop  a  procedure  which  will  differentially  augment  values  for  earthquakes  and  explosions 
with  the  larger  increase  for  earthquakes. 

There  was  an  m^  5.5  event  at  04:11:24.7  at  41.1  °N  143. 5°E  whieh  caused  problems  at  IASA, 
and  perhaps  at  NORSAR,  for  event  1  3.  The  Rayleigh  waves  from  this  other  event  arrived  at 
LASA  at  almost  the  same  time  as  those  from  event  13,  and  made  it  impossible  to  get  an  M  meas¬ 
urement.  At  NORSAR,  the  P  arrived  at  about  the  time  of  Love  waves  from  event  13.  The  P 
can  be  seen  by  checking  polarization  on  Fig.  Ill— 1 3 .  The  Rayleigh  wave  for  event  1  3  at  NORSAR 
looks  acceptable  but  it  is  at  the  same  time  as  PP  from  the  interfering  event.  This  eould  aeeount 
for  the  high  value,  but  is  not  likely. 

It  is  interesting  to  note  that  event  15  was  identified  as  problem  event  II  in  the  earlier  I  ASA 
report^  in  whieh  the  event  was  discriminated  by  IASA  Mg -m^,  although  it  did  have  a  somewhat 


small  Mg  value. 


R.  T.  Laeoss 


D.  INITIAL  DISCRIMINATION  RESULTS  AT  NORSAR 

An  initial  study  has  been  made  of  the  detection  and  classification  capabilities  of  NORSAR 
using  the  body-wave-to-surface-wave-magnitude  (m^-Mg)  criterion.  The  ehief  purpose  of  this 
study  is  to  identify  lines  of  future  research  in  discrimination  procedures  using  NORSAR  data. 

The  major  difficulty  in  applying  the  m^-Ms  criterion  is  the  measurement  of  surfaee-wave 
amplitudes  of  events  presumed  to  be  underground  explosions  —  the  presumption  being  based  on 
body-wave  detections  and  location.  By  using  an  array,  the  most  simple  way  to  increase  the  sig¬ 
nal  amplitude  with  respect  to  the  noise  field  is  to  align  the  traees  in  time  and  sum;  this  sum  or 
beam  can  then  be  used  for  the  measurement  of  surfaee-wave  amplitudes  if  the  noise  is  sufficiently 
suppressed.  In  this  study,  beamforming  is  the  main  technique  used  in  signal  enhancement.  In 
order  to  apply  this  technique,  especially  for  weak  signals  that  cannot  be  discerned  on  single 
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traces,  one  must  have  some  idea  of  the  rate  of  propagation  of  energy  from  the  source  region  to 
the  array  site  (i.e.,  the  group  velocity  along  a  given  path)  and  the  rate  of  propagation  of  energy 
across  the  array  (i.e.,  the  phase  velocity  at  the  array  site).  Additionally,  once  the  beam  is 
formed,  it  is  desirable  to  know  how  much  signal  loss  is  to  be  expected  from  the  beamforming 
process.  Because  of  the  dispersive  nature  of  surface  waves,  all  these  quantities  are  a  function 
of  period  and  propagation  path;  here,  I  have  concentrated  on  the  20-sec  components  from  sources 
in  central  Asia. 

Figure  III  —  1  5  shows  a  suite  of  Rayleigh-wave  group-velocity  measurements  made  at  NORSAR 
from  four  events  at  a  distance  of  about  40°  and  an  azimuth  of  about  75°.  A  clear  minimum  at 
about  the  18-sec  period,  typical  for  continental  regions,  is  evident  and  3.0  km/sec  has  been  taken 
as  a  characteristic  of  the  20-sec  group.  Phase-velocity  measurements  at  NORSAR  are  some¬ 
what  less  precise  because  of  the  limited  areal  extent  of  the  array  (100  km).  The  high-resolution 
frequency-wavenumber  technique  was  applied  to  two  events  from  central  Asia.  The  scatter  is 
quite  serious,  but  3.  5  km/sec  has  been  taken  as  an  appropriate  velocity  to  be  used  in  aligning  the 
20-sec  portions  of  the  Rayleigh-wave  recordings.  It  is  noted  here  that  the  direction  of  approach 
as  yielded  by  the  high-resolution  technique  was  always  within  10°  of  the  computed  direction  for 
these  sources  and  at  most  periods  between  10  and  50  sec,  within  5°. 

With  these  velocities  and  using  computed  great  circle  azimuths,  vve  can  proceed  to  investi¬ 
gate  the  effects  of  heamforming  on  the  Rayleigh-wave  recordings.  Figure  IIl-l6(a)  shows  the 
vertical  motion  of  the  Rayleigh  waves  from  a  presumed  explosion  in  Eastern  Kazakh  on  2  5  April 
1971.  The  traces  have  been  aligned  and  summed  and  the  beam  is  shown  as  the  topmost  trace. 

In  this  figure,  the  nature  of  the  group  velocity  curve  is  seen,  the  initial  long-period  energy  is 
followed  abruptly  by  the  shorter-period  components  which  are,  in  turn,  followed  by  20-sec 
arrivals.  Qualitatively,  as  seen  on  the  beam,  the  20-sec  arrivals  have  been  enhanced  relatively 
to  shorter  periods  due  to  the  beamforming  process.  Figure  III-l6(a)  quantifies  this  observation; 
here,  the  power  spectrum  of  the  beam  is  compared  with  the  average  spectrum  of  the  14C-ring 
long-period  vertical  sensors.  These  instruments  are  evenly  spaced  around  the  300-km  circum¬ 
ference  of  this  outer  ring.  The  spectra  are  computed  using  512  points  sampled  at  1-sec  intervals 
and  are  smoothed  by  taking  the  average  of  all  components  in  a  0.01-Hz  band.  In  the  right-hand 
side  of  Fig.  III-l6(b),  the  difference  between  the  beam  and  the  average  spectrum  is  plotted  and 
it  is  seen  that  in  this  single  case,  when  we  have  a  high  signal-to-noise  ratio,  the  beam  spectrum 
is  down  relative  to  the  array  by  about  2  dB  at  20  sec.  At  frequencies  higher  than  0.2  Hz,  the  beam 
spectrum  is  down  in  the  range  10  to  1  5  dB  below  the  array  average.  By  assuming  that  the  signal 
contains  no  energy  at  these  higher  frequencies,  this  can  be  considered  the  noise-suppression 
capability  of  the  beam  in  this  example.  In  any  case,  the  signal  deterioration  of  2  dB  (0.1  of  a 
magnitude  unit)  at  20  see  caused  by  beamforming  these  traces  is  the  relevant  factor  that  will  be 
applied  to  subsequent  beams. 

As  an  additional  example  of  the  beamforming  capability,  Fig.  III-17  shows  the  aligned  traces 
and  beam  for  another  presumed  explosion  from  the  same  region  as  the  example  of  Fig.  III-l6(a). 
Although  the  signal-to-noise  ratio  here  is  much  smaller  than  in  the  previous  example,  we  feel 
that  a  reliable  is  measurable  from  the  beam  even  though  the  arrival  is  not  discernible  on  the 
individual  traces. 
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A  preliminary  investigation  has  been  made  of  one  other  signal- enhancement  technique  during 
this  study.  Figure  III-18  shows  the  same  event  as  Fig.  IIl-l6(a).  However,  in  this  case  the  sec¬ 
ond  trace  from  the  top  has  been  replaced  with  a  180-sec  portion  of  the  top  trace  which  is  the  beam. 
The  second  trace  was  then  cross-correlated  with  the  top  trace,  and  the  third  trace  represents 
the  results  of  this  cross-correlation.  Although  180  sec  is  relatively  short,  an  experiment  has 
shown  that  doubling  the  length  of  the  reference  trace  (in  this  case  to  include  the  low-amplitude 
long-period  displacement)  increases  the  maximum  amplitude  of  the  cross -correlation  function 
by  only  10  percent.  This  same  reference  trace  of  Fig.  Ill— 1 8  was  applied  to  an  event  from  the 
same  region  but  in  a  case  where  1  would  decline  to  measure  an  Mg  for  the  array  beam.  This 
beam  is  shown  in  Fig.  111-19,  beneath  it  the  reference  trace  of  Fig.  Ill— 1 8 ,  and  beneath  that  the 
cross-correlation  function  of  the  two.  It  is  clear  from  the  latter  that  a  detection  has  been  made 
and  the  amplitude  can  be  measured  with  some  confidence.  The  left-hand  margins  of  Figs.  III-18 
and  -19  are  delayed  equally  with  respect  to  the  origin  times  of  the  individual  events. 

The  Mg-m^  criterion  has  been  applied  to  18  earthquakes  from  central  Asia  and  9  presumed 
Kazakh  explosions.  The  events  considered  are  at  distances  between  32°  and  64°  to  the  northeast 
of  NORSAR.  The  formula 

Mg  =  log10(A/T)  +  1.66  log  (A)  +  3.3  (III-l) 

was  applied  to  the  vertical  component  of  the  Rayleigh  wave,  where  A  is  the  maximum  amplitude 
in  microns  at  period  T,  and  A  is  the  distance  in  central  angle  degrees.  In  this  study,  A  was 
measured  at  periods  within  2  sec  of  20  sec.  With  the  exception  of  3  earthquakes,  all  the  m^  values 
were  taken  from  the  NOS’s  Preliminary  Determination  of  Epicenters  (PDE)  publication  for  events 
having  depths  less  than  38  km  or  reported  as  "normal"  depth.  In  Fig.  Ill -2 0,  the  results  of  plot¬ 
ting  Mg  vs  m^  are  shown.  The  magnitudes  of  the  earthquakes  numbered  8,  14  and  15  are  taken 
from  the  daily  bulletin  of  the  Seismic  Array  Analysis  Center  (SAAC).  In  Fig.  111-20,  we  have 
plotted  the  empirical  Mg-mb  relation  of  Gutenberg  and  Richter/3  which  is  Mg  =  1.59  —  3.97. 

Although  0.1  of  a  unit  of  surface-wave  magnitude  has  been  added  to  all  Mg  values  based  on 
beam  measurements,  this  correction  will  not  affect  the  separation  of  source  types  that  has  been 
achieved  in  Fig.  111-20.  The  Mg  of  explosion  point  9  was  measured  from  the  cross-correlation 
function  of  Fig.  Ill— 1  9.  The  most  surprising  aspect  of  the  figure  is  the  apparent  steepness  of  the 
fairly  well-defined  line  the  explosion  points  represent.  Equation  (III— 1 )  was  applied  to  earthquakes 
recorded  at  NORSAR  for  which  NOS  had  reported  an  Mg,  and  the  two  agreed  within  0.2  of  a  unit. 

It  is  difficult  to  attribute  the  slope  of  the  line  to  a  propagaxion-path  effect  since  all  these  explosion 
surface  waves  had  essentially  the  same  path  to  NORSAR.  .'herefore,  assuming  that  the  NOS 
determinations  arc  consistent  and  these  measurements  are  accurate,  the  steep  slope  of  the 
explosion  line  must  be  attributed  to  a  source  effect.  The  implication  is  that  the  short-  and  long- 
period  displacements  at  the  source  do  not  scale  equally  with  yield. 

Explosion  data  from  the  Ural  Mountains  region  have  not  been  discussed  here  since  they  fall 
within  a  range  where  the  direct  application  of  Eq.  (III-l)  is  not  valid,  and  at  present  there  is  not 
a  significant  population  of  earthquakes  against  which  to  test.  This  should  be  a  line  of  further 
research.  It  is  impossible  to  make  any  general  valid  statements  about  the  capability  of  an  array 
after  such  a  short  study.  However,  one  point  becomes  evident;  using  elementary  techniques  of 
data  processing  described  here,  the  lower  threshold  of  surface-wave  detection  has  not  been 
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reached  when  working  with  events  near  the  lower  threshold  of  NOS-reported  location  capability 
in  central  Asia.  In  order  to  establish  the  ultimate  long-period  capability  of  NORSAR,  other 
reliable  means  of  short-period  detection  and  location  of  events  in  central  Asia  must  be  used  or 

devel°Ped-  J.  R.  Filson 

E.  LASA  HIGH-VELOCITY  CORRECTIONS 

One  problem  that  has  always  plagued  beamforming  with  IASA  is  the  large  time  anomalies 
that  exist  within  the  array.  Without  sufficient  knowledge  of  these  corrections,  it  would  not  be 
possible  to  obtain  correct  or  maximum  beams.  The  time  corrections  for  the  short-period  P-phase 
were  determined  some  time  ago  and  have  been  in  use  for  several  years.  Only  recently  has  em¬ 
phasis  been  placed  on  beamforming  for  higher -velocity  phases  sueh  as  PKP  or  core  phases  such 
as  PKiKP  or  PKJKP.  The  earlier  analysis  of  the  LASA  time  anomalies  indicated  that  the  core¬ 
phase  corrections  were  large,  often  larger  than  the  P-wave  anomalies.  With  the  aid  of  data  ob¬ 
tained  from  IBM-SAAC  and  NOs/  these  high-veloeity  eorreetions  have  been  determined  and  are 
now  in  use  in  our  PDP-7  beamforming  programs. 

The  initial  data  contained  about  200  core  phases  that  were  processed  by  the  event  processor 

at  SAAC.  One  of  the  outputs  of  the  event  processor  is  a  list  of  measured  delay  times  for  the  phase 

as  seen  at  IASA.  These  measured  times  were  compared  with  the  computed  plane-wave  times  for 

the  phase  to  yield  the  plane-wave  anomalies  whieh  were  tabulated  into  cells  that  cover  each  10° 

in  azimuth  and  0.02  sec/deg  in  slowness.  In  the  range  of  interest,  0  to  4see/deg,  most  cells 

did  not  contain  any  data  points.  It  was  necessary  to  approximate  these  data  with  a  function  which 

could  then  be  evaluated  at  each  cell.  Earlier  experience  with  the  P-wave  corrections  indicated 

£ 

that  the  time  anomalies  could  be  approximated  by  a  two-term  Fourier  series  with  azimuth.  Fur¬ 
ther,  the  core  anomalies  seemed  to  have  a  second-  or  third-order  polynomial  dependence  with 
distance  (slowness). 

Defining  the  station  correction  as  a  two-term  Fourier  series  in  azimuth  and  a  third-order 
polynomial  in  slowness  meant  that  some  twenty  coefficients  were  neeessary  to  express  the  func¬ 
tion.  However,  when  the  condition  is  imposed  that  requires  the  function  to  be  a  constant  (inde¬ 
pendent  of  azimuth)  for  zero  slowness,  the  number  of  coefficients  reduces  to  eleven.  All  data 
for  each  subarray  have  been  least-square  fitted  to  the  restrained  function  defined  by  using  the 
eleven  coefficients.  The  resulting  equation  was  evaluated  at  each  cell  point  and  compared  re¬ 
markably  well  with  the  original  data.  Each  cell  point  now  had  a  data  value  obtained  either  from 
real  data  or  from  the  equation.  In  order  to  obtain  a  more  uniform  table  of  values,  the  table  was 
smoothed,  requiring  that  the  real  data  points  would  remain  unchanged  while  the  function  values 
eould  change.  The  resulting  tables  of  high-velocity  LASA  station  corrections  arc  now  in  use  in 
all  beamforming  programs  and  are  yielding  excellent  results,  as  indicated  in  the  successful 
search  for  the  eore-phase  PKJKP  described  in  Sec.  Vll. 

R.  M.  Sheppard 
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Fig.  111-1.  Beam  pattern  for  long-period  array  at  NORSAR. 
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Fig.  1 1 1  —2  (a  -c) .  High-resolution  frequency-wavenumber  spectra  for  LPZ  noise  sample  at  NORSAR. 
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Fig.  Ill-3(a-c).  High-resolution  frequency-wavenumber  spectra  for  LPNS  noise  sample  at  NORSAR. 
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Fig.  Ill-4(a-c).  High-resolution  frequency-wavenumber  spectra  for  LPEW  noise  sample  at  NORSAR. 
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Fig.  1 11-5.  Power  spectra  for  long-period  noise  sample  at  NORSAR. 
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Fig.  MI-6.  Beam  pattern  far  ALPA. 
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Fig.  111-7.  Location  of  16  September  1970  noise  source  at 


0640  to  0740  using  ALPA  and  LASA  simultaneously. 
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Fig.  Ill -8.  Locotion  of  20  October  1970  noise  source  ot  0550  to  0650  using  ALPA  ond  LASA  simultaneously. 
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Fig.  IM-9.  Comparison  of  power  spectra  of  noise  at  (a)  ALPA  and  (b)  LASA  from  0640  to  0740 
on  16  September  1970.  Power  is  measured  in  decibels  relative  ta  1  mp^/Hz  at  0.04  Hz. 
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Fig.  111-11.  NORSAR  configuration  with  outline  of  geology. 
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Fig.  111-12.  NORSAR  long-period  vertical  records  for  four  Asian  events.  All  signals  normalized  to  peak  values. 


Section  111 


39 


Fig.  111-13.  Three-component  long-period  records  for  fourteen  Asian  events  recorded  in  Norway. 
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Fig.  111-14.  Ms-rr»b  plot*  using  NOS  (o)  Norway  Rayleigh-wave  M$ with  period 
in  18-  to  23-sec  range;  (b)  Norway  Love-wave  M$  far  periods  greater  than  18  sec; 
(c)  larger  value  from  (o)  ond  (b)  used  for  M$;  and  (d)  larger  value  fram  (c)  and  LASA 
used  far  M$. 


Fig.  111-15.  Rayleigh-wave  graup-velacity 
measurements  ot  NORSAR  fram  4  events  in 
central  Asia. 
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Fig.  1 1 1  —  1 6  (a) .  Beam  and  aligned  traces  of  Rayleigh  waves  fram  Eastern  Kazakh  event 
af  25  April  1971. 
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Fig.  1 1 1  - 1 6 (b) .  Power  spectrum  of  beam  and  array  average-pawer  spectrum,  and  difference 
between  the  twa.  Power  is  measured  relative  to  1  m^i  2/ Hz  at  20  sec. 
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Fig.  111-17.  Beam  and  aligned  traces  af  Eastern  Kazakh  event  af  25  May  1971. 
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Fig.  111-18.  Rayleigh-wave  beam  af  event  of  25  April  1971;  reference  waveform  fram  90-sec 
portion  af  this  beam;  and  cross-correlatian  of  reference  trace  and  beam. 
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Fig.  111-19.  Rayleigh-wave  beam  af  event  fram  Eastern  Kazakh  an  9  October  1971;  reference 
waveform  fram  event  of  Fig.  111-18;  and  cross-correlation  af  reference  trace  and  beam. 


43 


Section  III 


Fig.  111-20.  Ms-m^  diagram.  M$  measured  at  NORSAR,  m^  NOS  or  SAAC  determination. 
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IV.  NETWORK  CHARACTERISTICS 


A.  DETECTION  AND  DISCRIMINATION  OF  CENTRAL  ASIAN  EARTHQUAKES 
AND  EXPLOSIONS  WITH  A  CONTINENTAL-SIZE  WWSSN  NETWORK 

A  study  has  been  undertaken  to  determine  the  characteristics  of  long-period  radiation  from 

1-3 

Central  Asian  earthquakes  and  explosions.  Such  a  study  will  supplement  previous  work  with 
digital  data  from  NORSAR,  ALPA  and  LPE.  At  present,  the  data  are  based  on  the  visual  analy¬ 
sis  of  World-Wide  Standard  Seismograph  Network  (WWSSN)  film  chips  contained  in  the  Special 
Event  Library  at  Lincoln  Laboratory.  In  addition  to  determining  the  detection  and  discrimination 
capability  of  an  Eurasian  WWSSN  network,  a  further  goal  of  this  study  is  to  anticipate  methods 
useful  in  optimizing  large-array  data  analyses  for  events  in  Central  Asia. 

Here  the  results  to  date  are  reported  for  data  collected  from  twelve  Asian  and  European 
WWSSN  stations.  Each  station  was  selected  on  the  basis  of  high  gain,  low  noise  and/or  strategic 
location.  A  list  of  the  stations  is  given  in  Table  IV- 1.  The  current  data  base  is  comprised  of 
21  presumed  Asian  nuclear  explosions  and  17  Central  Asian  earthquakes. 

A  test  of  the  standard  M  -m,  discriminant  was  made  from  the  events  listed  in  Table  IV- 2. 

s  b 

Figure  TV-1  shows  the  results.  The  body-wave  magnitude  was  that  given  by  the  National 
Ocean  Survey  (NOS)  when  three  or  more  stations  reported.  Otherwise,  the  NOS  value  was  com¬ 
bined  with  a  value  determined  from  our  12-element  network.  Even  when  a  sufficient  number  of 
NOS  stations  reported,  the  body-wave  magnitude  from  our  near-field  network  was  invariably 


TABLE  IV-1 

STATIONS  USED  IN  CENTRAL  ASIAN  STUDY 

Station  Symbol 

Station 

Latitude 

Longitude 

CHG 

Chiengmai,  Thailand 

1  8°47,24,IN 

98°58 ' 37"E 

KBL 

Kabul,  Afghanistan 

35°32'27"N 

69°02 ' 35"E 

KEV 

Keva,  Finland 

69°45 ' 00 1 '  N 

27°00'00"E 

LAH  (NIL,  1970) 

Lahore,  Pakistan 

31  °33‘  00"N 

74°20'00"E 

MSH 

Mashed,  Iran 

36°1 8'40"N 

59°35 1 1 6"E 

ND! 

New  Delhi,  India 

28°42'00,,N 

77°13'00"E 

NUR 

Nurmijarvi,  Finland 

60°30'30,,N 

24°39'07"E 

SHI 

Shiraz,  Iran 

24°38,40"N 

52°31 '34"E 

SHL 

Shillang,  India 

25°34'00MN 

91  °53'00"E 

POO 

Paana,  India 

1 8°32 ' 00" N 

73°51 '00"E 

QUE 

Quetta,  Pakistan 

30°1 1  1 03nN 

66°57 1 00"E 

TAB 

Tabriz,  Persia 

38°04,03"N 

46°19'36"E 

45 
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within  one- half  a  magnitude  unit.  The  surface- wave  magnitude  M  was  computed  by  measuring 

the  peak-to-peak  amplitude  A  (in  microns)  and  the  corresponding  period  T  (in  seconds)  of  the 

4 

Rayleigh  wave  at  its  maximum  and  applying  : 

Mg  =  log  ( A/T)  +  0.92  +  log  (A)  A  <  25° 

Mg  =  log  (A/T)  +  1.66  log  (A)  A  ^25° 

Maximum  values  were  read  to  maintain  maximum  detectability. 

For  a  given  event,  the  plotted  Mg  value  was  averaged  over  the  values  determined  at  the 
various  detecting  stations.  Typical  single-event  scatter  of  Mg  is  presented  in  Fig.  IV- 2.  When 
a  surface  wave  was  not  detected,  the  amplitude  of  the  noise  in  the  signal  band  was  measured. 

In  these  cases,  an  upper  limit  on  the  surface-wave  magnitude  was  computed  by  assuming  de¬ 
tectability  of  a  waveform  when  the  signal-to-noise  ratio  was  one-to-one.  Data  points  computed 
in  such  a  manner  are  marked  with  an  arrow.  The  tail  of  the  arrow  marks  the  average  of  the 

four  lowest  values,  while  the  head  of  the  arrow  points  to  the  single  lowest  value. 

Observations  based  on  the  M  -m,  measurements  are: 

s  b 

(1)  The  surface  waves  from  all  events  with  Ms  ^3.5  are  detected.  It  seems 
probable  that  the  inclusion  of  more  data  will  lower  that  figure  to  Ms  >  3.1. 

Further,  a  factor-of-four  improvement  in  signal  detectability  would  lower 
the  value  to  2.9  and  possibly  down  to  2.4.  It  seems  reasonable  that  digi¬ 
tization  of  the  network  would  give  such  an  improvement. 

(2)  The  events  clearly  separate  into  two  populations.  However,  any  linear 
relationship  between  Ms  and  m^  that  extends  below  m^  =5.0  will  have 
large  scatter.  The  most  troublesome  events  at  =  5.0,  Mg  =3.5  are 
geographically  confined  to  a  small  area  in  Tibet. 

The  separation  between  natural  and  artificial  M  -m,  values  has  been  shown  to  improve  by 

2  3  s  D 

measuring  magnitudes  at  longer  periods.  '  Since  we  used  maximum  Rayleigh-wave  amplitudes, 
the  poor  separation  at  low  magnitudes  shown  in  Fig.  IV- 1  might  be  improvable.  Two  of  the 
natural  events  in  our  data  base  contained  sufficient  broadband  energy  to  make  an  estimate  of 
the  decay  of  Mg  with  increasing  period.  A  plot  of  the  results  in  Fig.  IV- 3  indicates  that  a  drop 

of  at  least  one  unit  of  surface- wave  magnitude  can  be  expected  between  12  and  30  sec.  Such  a 

2 

decrease  is  consistent  with  North  American  explosion  data.  Even  though  our  data  are  quite 
limited,  the  supposition  of  a  much  improved  Mg-m^  discriminant,  based  on  longer  period  meas¬ 
urements,  seems  tenuous. 

Of  particular  interest  to  NORSAR  capabilities  are  the  characteristics  of  waves  observed  at 
the  Finland  stations,  Kevo  and  Nurmi jarvi.  As  shown  in  Fig.  IV- 2,  at  least  one  of  the  stations 
has  the  maximum  Mg  recorded  by  the  array  in  6  out  of  the  8  events,  often  by  half  a  unit  of 
surface- wave  magnitude.  In  the  15  events  in  the  complete  data  base  in  which  at  least  one  of  the 
stations  made  an  observation,  one  or  both  was  a  maximum  or  well  above  average  13  times.  The 
average  period  of  the  arrivals  was  12  to  14  sec.  It  seems  clear  that  the  continental  path  from 
Mid-Asia  to  Scandinavia  is  an  efficient  one  for  short-period  surface  waves  over  teleseismic  dis¬ 
tances.  Since  it  is  not  uncommon  to  have  relatively  low  noise  in  the  10-  to  15-sec  band,  NORSAR 
surface-wave  detection  capabilities  may  be  aided  by  examining  data  in  that  period  range.  It  has 

been  noticed  that  NORSAR  M  values  based  on  20-sec  Rayleigh-wave  data  from  Central  Asia  are 
5  S 

anomalously  low. 
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In  the  future,  NORSAR,  ALPA  and  LPE  data  will  be  included  in  our  analysis  and  the  com¬ 
bined  capabilities  of  these  arrays  with  WWSSN  will  be  used  in  determining  the  source  and  path 

characteristics  for  Eurasian  events.  m  „  T  , 

T.  E.  lenders 

B.  SHORT-PERIOD  NETWORK  CAPABILITY 

In  conjunction  with  the  experiment  mentioned  in  the  previous  section,  we  are  investigating 
issues  connected  with  the  production  of  a  more  adequate  low-magnitude  bulletin.  One  such  issue 
is  the  detection  level  and  credibility  of  the  Ij^SA  detection  processor;  another  is  the  detection 
level  of  the  WWSSN  stations  for  events  in  the  Sino- Soviet  area. 

The  experiment  is  being  done  in  two  phases.  The  first  phase  is  to  determine  the  world 
seismic  activity  for  a  given  period  -  in  this  case,  January  1970.  The  second  phase  will  then  be 
to  see  if  the  long-period  waves  from  the  events  can  be  detected  on  the  worldwide  network.  An 
important  aspect  of  this  experiment  is  to  set  it  up  in  such  a  way  that  the  procedures  used  may 
be  suitable  for  the  routine  processing  of  any  or  all  periods  in  the  future. 

In  an  effort  to  compile  as  complete  a  list  of  epicenters  as  possible  for  this  study,  the  National 
Ocean  Survey  recomputed  all  their  data  for  January  19  70.  The  more  than  700  events  they  were 
able  to  locate  will  form  one  of  the  main  data  sources.  For  this  same  month,  the  LASA  detection 
log  indicated  more  than  2000  possible  event  detections.  The  input  from  these  two  sets  of  data 
will  form  the  main  starting  point  for  the  effort  to  produce  as  complete  a  list  of  epicenters  as 
possible.  Since  many  of  the  events  in  the  initial  list  will  be  phases  or  false  alarms,  it  will  be 
necessary  to  investigate  each  entry  individually.  Those  events  that  are  phases  can  usually  be 
eliminated  in  the  first  step.  When  this  is  done,  the  NOS  list  will  yield  epicenters  with  a  high 
degree  of  reliability.  However,  the  IASA  detection  log  is  left  with  a  large  number  of  detections 
that  must  be  individually  investigated.  An  example  of  this  is  an  event  which  appeared  on  the 
LASA  bulletin  as  an  unlocated  event  arriving  from  the  east.  This  event  recorded  best  on  a  beam 
centered  at  37"  N  26°  E  and  triggered  an  additional  38  beams.  By  using  the  location  of  the  beam 
as  the  epicenter,  arrival  times  were  computed  to  some  of  the  near  WWSSN  stations.  This  event 
could  later  be  seen  on  ATU  (Athens  University,  Greece)  and  EIL  (Eilat,  Israel)  as  shown  in 
Fig.  IV- 4.  The  epicenter  was  located  at  37.  2°  N  28.1°  E,  and  had  a  magnitude  less  than  4.0. 

Each  LASA  detection  log  event  is  being  studied  in  this  way.  Early  results  indicate  we  can  expect 
that  50  or  60  events  per  day  will  be  verified  in  this  way.  The  addition  of  a  detection  log  from 
any  other  large  array,  such  as  NORSAR,  may  double  this  number. 

R.  E.  Needham 
R.  M.  Sheppard 
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Fig.  IV—  1  -  Ms-m[3  relationship  far  Asian  earthquakes  and  explosions.  Arrows  indicate 
that  an  upper-limit  estimate  based  an  signal  band  noise  was  made.  Maximum  Rayleigh- 
wave  amplitudes  were  used.  Letter  symbols  identify  events;  number  af  recording 
stations  far  each  event  is  indicated  by  numerals  in  parentheses. 
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Fig.  IV-2.  Single-event  Ms  scatter.  Event  symbols  refer  to  Table  IV-2.  Note  behavior  of  Finland  Stations,  Nurmijarvi  and  Kevo. 


Section  IV 


Fig.  IV-3.  Behavior  af  M$  as  a  function  af  period  far  two  Asian  earthquakes.  Event  symbols 
refer  ta  Table  IV-2. 
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Fig.  IV-4.  Verification  of  a  LASA  beam  detection  on  shart-period  seismometers 
at  Athens  and  Eilat. 
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V.  LONG-PERIOD  SEISMOMETERS 


A.  NONPROPAGATING  NOISE  IN  LONG-PERIOD  SEISMOMETERS 

It  has  been  known  for  some  time  that  the  noise  measured  by  long-period  seismometers  is 

composed  of  two  parts:  propagating  noise  which  is  coherent  over  distances  of  at  least  100  km, 

and  nonpropagating  noise  which  is  incoherent  over  distances  greater  than  a  few  kilometers.  The 

nonpropagating  component  is  therefore  completely  uncorrelated  between  long-period  instruments 

at  LASA  and  can  only  be  reduced  by  by  summing  n  instruments.  However,  Capon"*  has  shown 

using  coherence  spectra  that  there  is  a  significant  correlation  between  measured  long-period 

seismic  noise  and  microbarograph  recordings  at  the  same  site  at  times  when  the  nonpropagating 

component  is  appreciable.  Since  this  component  contributes  at  least  40  percent  of  the  noise  up 

1 

to  50  percent  of  the  time,  it  seemed  worthwhile  to  consider  the  possibility  of  eliminating  this 
noise  using  the  fact  that  it  is  correlated  with  the  output  of  the  nearby  microbarograph. 

We  can  divide  the  seismometer  output  S(t)  into  two  parts:  that  which  is  correlated  with  the 
microbarograph  output,  and  therefore  caused  by  local  atmospheric  pressure  fluctuations  N{t); 
and  that  part  not  correlated  with  the  microbarograph  output  E(t).  Thus, 

S(t)  =  N(t)  +  E(t)  .  (V-l) 

It  is  possible  to  define  an  atmospheric  pressure  fluctuation  in  the  vicinity  of  the  seismometer 
and  microbarograph  which  would  produce  a  unit  positive  impulse  in  the  output  of  the  micro¬ 
barograph,  and  some  other  response  I(t)  in  the  output  of  the  seismometer.  In  this  case,  if  the 
microbarograph  output  is  M{t),  the  pressure-induced  noise  N(t)  in  the  seismometer  output  is 
given  by 

N(t)  -  C  M(T)I(t-T)dr  .  (V-2) 


If  we  define  the  cross-correlation  function  to  be 


r 

c  —  o 


^ABU)  =  '  A(t  +  n  B(t)  dt 


it  can  be  shown  that 


■r 

-0 


1(3  +  T)  %M(S)  ds 


( V-3) 


( V-4) 


where  wc  have  used  the  knowledge  that  E{t)  and  M(t)  are  uncorrelated  (by  definition)  or, 
mathematically, 

<^EM(r)  =  0  for  all  t  .  (V-5) 

Equation  (V-4)  can  be  solved  for  I(t),  which  can  be  convolved  with  the  microbarograph  output  to 
obtain  an  estimate  of  the  pressure-induced  noise  [Eq.  ( V-2 )  1 .  This  predicted  noise  can  be  sub¬ 
tracted  from  the  seismic  signal  S(t)  to  reduce  the  noise  level. 

This  method  was  applied  to  a  noise  sample  from  the  3-componcnt  instrument  at  A0  at  LASA, 

1 

beginning  17h20m00s  on  23  August  1967.  Capon  showed  that  the  nonpropagating  component 
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contributed  nearly  10  dB  more  to  the  noise  measured  by  the  vertical  component  than  did  the 
propagating  part.  Figure  V-l  shows  the  microbarograph  output  (1),  the  seismometer  output  (2) 
the  predicted  noise  (3),  and  the  seismometer  output  with  the  predicted  noise  subtracted  (4)  for 
each  of  the  three  components  (a)  vertical,  (b)  north-south,  and  (c)  east-west.  In  Fig.  V-2  can 
be  seen  the  coherence  spectra  between  the  microbarograph  and  the  seismometer  (1),  and  be¬ 
tween  the  microbarograph  and  the  seismometer  output  with  the  predicted  noise  subtracted  (2) 
for  (a)  vertical,  (b)  north-south  and  (c)  east-west  components.  The  big  reduction  in  coherence 
in  the  seismic  signal  bandwidth  from  Fig.  V-2(l)  to  2(2)  for  each  component  indicates  that  much 
of  the  pressure-induced  noise  has  been  eliminated.  It  is  therefore  possible  to  reduce  this  kind 
of  noise  significantly  at  every  site  at  which  there  is  also  a  microbarograph. 

This  work  does  not  indicate  how  much  of  the  pressure-induced  noise  is  due  to  buoyancy  ef¬ 
fects  in  the  seismometer,  or  how  much  is  the  result  of  ground  displacements  caused  by  fluctua¬ 
tions  in  atmospheric  pressure.  The  atmospherically  produced  ground  displacement  is  a  spatial 

integration  of  the  whole  pressure  pattern  around  the  seismometer;  the  fluctuations  in  ground 

2 

displacement  depend  on  the  variation  of  the  pressure  pattern  with  time.  We  would  therefore 
expect  the  stability  of  the  function  I(t)  to  depend  on  the  weather  if  the  ground-loading  effect  were 
dominant.  On  the  other  hand,  if  the  buoyancy  effect  were  dominant,  we  would  expect  I(t)  to  be 
as  stable  as  the  responses  of  the  microbarograph  and  seismometer. 

To  determine  which  of  these  effects  is  the  more  important,  the  previously  computed  transfer 
functions  I(t)  for  each  seismometer  component  at  AO  were  applied  to  a  noise  sample  beginning 
19h50m00s  on  1  July  1967,  which  Capon  found  to  have  a  large  nonpropagating  component. 

The  microbarograph  output  (1),  the  seismometer  output  (2),  the  predicted  noise  [using  the  I(t) 
calculated  from  the  23  August  sample]  (3),  and  the  difference  between  the  seismometer  output 
and  the  predicted  noise  (4)  are  shown  for  each  component  of  the  seismometer  at  AO  in  Fig.  V-3. 
The  predicted  noise  appears  to  bear  little  relation  to  the  measured  noise,  and  on  subtraction 
from  it  even  increases  the  noise  amplitude.  We  may  conclude  from  this  and  the  foregoing  argu¬ 
ment  that  the  ground-loading  effect  is  dominant. 

3 

Herrin,  et  aJ.  have  reported  coherencies  between  noise  measured  by  long-period  seis¬ 
mometers  and  microbarographs  at  Grand  Saline,  Texas.  They  were  particularly  careful  to  iso¬ 
late  the  seismometers  from  buoyancy  effects,  and  their  results  clearly  show  that  the  measured 

2 

noise  exhibited  the  qualities  which  would  be  expected  from  Sneddon’s  theoretical  considerations  : 
the  decay  of  amplitude  with  depth  of  seismometer  emplacement,  and  the  dependence  on  the  way 
in  which  the  surrounding  pressure  pattern  varies. 

Since  the  ground-loading  effect  has  a  critical  spatial  as  well  as  time  dependence,  the  noise 
measured  by  a  seismometer  could  be  predicted  with  some  confidence  only  from  an  array  of 
microbarographs  all  situated  within  a  radius  of  a  few  kilometers  of  the  seismometer. 

A.  M.  Ziolkowski 

B.  DISCUSSION  OF  DATA  FROM  THE  IGPP  QUARTZ  ACCELEROMETER 

4 

Long-period  digital  data  from  a  new  wideband  accelerometer  have  been  supplied  to  us  by 
Professor  B.  Block  of  IGPP,  University  of  California,  San  Diego.  The  instrument  was  located 
at  a  surface  site  (vault  at  a  depth  of  10  feet)  near  the  coast  on  the  boundary  of  a  Naval  Air  Base 


56 


Section  V 


in  California  (32.88°N,  117.1^).  The  output  of  primary  interest  to  us  had  a  response  curve 
roughly  matched  to  that  of  instruments  used  by  Pomeroy,  et  al.5  The  shape  of  the  response  is 
shown,  along  with  the  LP  LASA  response  and  the  WWSSN  response,  in  Fig.  V-4.  Henceforth, 
reference  to  the  IGPP  instrument  will  mean  the  output  achieved  through  the  filters  which  achieve 
this  response  curve. 

Figure  V-5(b)  shows  the  results  of  noise  spectra  obtained  from  the  IGPP  instrument  during 
April  and  May  1971.  No  correction  for  instrument  response  has  been  made.  About  90  percent 
of  the  observations  lie  within  the  dashed  lines.  We  would  guess  that  the  height  of  the  peaks  at 
0.06  and  0.1 2  Hz  would  be  larger  and  the  entire  noise  level  raised  during  late  fall  and  winter 
months,  but  such  data  are  not  available.  Figure  V-5(a)  shows  the  spectra  of  Rayleigh  waves 
for  two  events  that  are  typical  of  many  which  have  been  measured.  From  these  spectra,  we 
note  that  the  best  signal-to-noise  ratio  exists  in  the  0.03-  to  0.05-Hz  band.  This  agrees  with 
results  obtained  using  several  narrowband  filters  at  IGPP^  and  at  Lincoln  Laboratory. 

We  recognize  that  sites  with  less  noise  at  frequencies  in  the  range  0.03  to  0.02  Hz  can  be 

found,  which  will  increase  the  signal-to-noise  ratio  at  those  longer  periods.  However,  data  in 
7 

the  literature  indicate  to  us  that  detection  capability  for  earthquakes  will  not  be  better  at  the 
longer  periods  at  those  sites  and  it  will  be  worse  for  explosions.  Also,  good  continental  paths 
can  exhibit  a  strong  phase  with  periods  as  short  as  10  sec  (see  Sec.  IV-A).  Such  a  phase  may 
occasionally  achieve  higher  signal-to-noise  ratios  than  in  the  0.03  to  0.02  band,  but  that  should 
normally  be  the  exception. 

Instruments  with  capabilities  similar  to  the  IGPP  instrument  have  often  been  favorably 
compared  with  WWSSN  instruments  for  detection  purposes,  and  it  is  often  assumed  this  is  be- 

7 

cause  the  newer  instruments  are  not  subject  to  atmospheric  buoyancy  near  40-sec  periods. 

This  is  probably  partly  correct  but  it  is  interesting  to  consider  what  would  happen  if  the  IGPP 
instrument  output  were  viewed  through  a  noiseless  filter  designed  to  give  the  same  frequency 
characteristics  as  a  WWSSN  station.  From  Fig.  V-4  we  see  that  at  8  sec  the  IGPP  instrument 
is  down  almost  30  dB  relative  to  the  WWSSN.  Assuming  the  noise  spectra  of  Fig.  V-5(b)  rep¬ 
resent  ground  motion  and  not  buoyancy,  we  conclude  that  if  the  IGPP  instrument  were  corrected 
to  the  WWSSN  response,  the  noise  at  8  sec  would  be  about  20  dB  above  the  noise  at  40  sec,  and 
10  to  15dB  above  the  noise  at  20  to  25  sec.  A  signal  received  under  those  conditions  would  re¬ 
quire  a  very  strong  signal-to-noise  ratio  in  the  20-  to  40-sec  band  to  be  detected  on  a  trace 
adjusted  so  that  the  8-sec  microseisms  would  not  be  too  large. 

Figure  V-6(a-b)  graphically  demonstrates  the  situation.  A  signal  from  Greece  was  re¬ 
corded  by  the  IGPP  instrument  with  a  good  signal-to-noise  ratio.  A  high-pass  filter  was  applied 
to  reduce  the  25-sec  signal  by  about  23  dB  and  not  affect  the  8-sec  noise.  Roughly,  this  gives 
a  trace  with  about  the  same  relative  amounts  of  power  at  8  and  25  sec  which  would  be  seen  with 
a  WWSSN  response  curve.  The  filtered  signal  is  barely  visible  and  we  can  see  why  a  WrWSSN 
record  would  show  poor  detection  capability  even  if  the  instruments  were  not  influenced  by 
buoyancy. 

This  suggests  that  WWSSN  station  detection  capability  could  be  significantly  enhanced  by 
using  digital  recordings  with  large  dynamic  range  to  allow  subsequent  filtering,  or  by  directly 
generating  filtered  records  with  the  8-  and  16-sec  microseisms  considerably  attenuated.  Of 
course,  if  buoyancy  effects  at  20-  to  30-sec  periods  are  too  large,  the  gain  would  not  be  as 
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large  as  indicated  by  the  above  demonstration.  We  believe  that  relatively  modest  isolation 
from  atmospheric  pressure  would  be  sufficient  most  of  the  time.  This  statement  is  based  upon 
a  consideration  of  LASA  data  obtained  from  instruments  with  limited  pressure  isolation.  Al¬ 
though  there  apparently  are  days  when  LASA  may  be  significantly  influenced  by  buoyancy  in  the 
20-  to  40-sec  range,  it  typically  constitutes  no  more  than  about  half  of  the  recorded  amplitude, 
and  even  some  of  that  must  be  due  to  tilt  which  is  ground  motion  and  cannot  be  removed  in  any 
simple  way  (see  Sec.  A  above,  and  also  Ref.  1).  Also,  atmospheric  problems  are  more  severe 
in  the  30-  to  40- sec  band  than  in  the  20-  to  30-sec  band.  Fall  and  winter  days  with  large  micro¬ 
seisms  would  probably  be  even  less  favorable  to  the  WWSSN  recordings. 

We  cannot  fairly  discuss  the  situation  which  would  obtain  at  a  site  very  quiet  at  40  sec  where 
the  buoyancy  problem  would  be  more  severe  except  for  instruments  such  as  the  IGPP  accelerom¬ 
eter  or  standard  instruments  extremely  isolated  from  the  environment,  such  as  the  Lamont  LPE 

installations.  0  m  T 

R.  1 .  Lacoss 

E.  N.  Ashley,  Jr. 
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Fig.  V-1.  Data  from  LASA  site  AO  beginning 
17h20m00s  an  23  August  1967.  Far  each  af 
three  components  of  seismometer  (a)  vertical, 
(b)  narth-south,  and  (c)  east-west  are  shown 
seismometer  output  (2),  predicted  noise  in 
that  output  (3),  and  seismometer  output  minus 
predicted  naise  (4).  In  each  case,  (1)  is 
microbarograph  output  at  AO. 
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Fig.  V-2.  From  data  shown  in  Fig.  V-l  are  the  computed  coherencies  between  microbarograph  and  seismometer  outputs  (1),  and  microbarograph 
output  and  seismometer  output  with  predicted  noise  subtracted  (2)  for  each  component  of  seismometer  (a)  vertical,  (b)  north-south,  and  (c)  east-west. 
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(Q)  VERTiCAL  COMPONENT 


Fig.  V-3.  Doto  from  LASA  site  AO  beginning 
19h50m00s  on  1  July  1967.  As  in  Fig.  V-l 
for  eoch  of  three  components  of  seismometer 
(o)  verticol,  (b)  north-south,  ond  (c)  eost- 
west  ore  shown  seismometer  output  (2),  pre¬ 
dicted  noise  in  thot  output  (3),  ond  seismom¬ 
eter  output  minus  predicted  noise  (4).  In 
eoch  cose,  (1)  is  microborogroph  output,  but 
predicted  noise  hos  been  computed  from 
microborogroph  output  using  transfer  functions 
colculoted  from  dota  of  Fig.  V-l* 
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(C)  EAST-WEST  COMPONENT 
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Fig.  V-4.  Relative  response  curves  of  long-period  instruments  of  WWSSN,  LASA 
and  (GPP  approximation  to  Lamont  LPE, 


62 


Section  V 


63 


Section  V 


VI.  DATA-PROCESSING  FACILITIES 


A.  LINCOLN  COMPUTER  FACILITIES 

Floating-point  hardware  has  been  ordered  for  the  two  PDP-7  computers.  When  installed, 
this  hardware  will  speed  up  by  an  order  of  magnitude  floating-point  operations  (now  done  by  soft¬ 
ware  routines  using  the  fixed-point  hardware).  The  Fortran  system  is  being  reprogrammed  to 
take  advantage  of  this  hardware.  The  new  internal  arithmetic  library  subroutines  will  operate 
about  5  times  faster,  and  the  science  library  (sines,  exponents,  etc.)  will  operate  about  50  times 
faster.  Machine  language  programs  can  take  full  advantage  of  this  hardware  and  speed  up  a 
floating  operation  by  about  10  times  over  the  old  Fortran  system. 

A  half-hour  16-mm  sound  film  has  been  made  of  the  Data  Analysis  Console  demonstrating 
its  various  features.  This  film  is  available  to  anyone  upon  request.  A  completely  revised  report 
on  the  Analysis  Console  is  also  being  prepared  and  will  be  available  to  those  interested. 

P.  L.  Fleck,  Jr. 

B.  IMPROVED  LIBRARY  ACCESS 

Because  of  the  large  numbers  both  of  magnetic  tapes  (about  10,000)  that  have  been  acquired 
over  the  last  few  years,  and  of  epicenters  that  are  reported  annually  by  NOS,  SAAC  and  the  Inter¬ 
national  Seismological  Center,  the  process  of  associating  epicenter  data  with  our  tape  library 
is  being  adapted  to  the  Lincoln  Information  Storage  and  Associative  Retrieval  System  (LISTAR).^ 
This  is  an  on-line  system  which  will  permit  the  user  to  define,  search,  modify  and  cross¬ 
associate  data  files.  In  our  case,  the  two  data  files  will  contain  epicenter  data  from  all  sources 
and  magnetic-tape  data  such  as  the  tape  start  time,  array  name  and  tape  location.  Input  of  seis¬ 
mic  data  into  this  system  will  provide  a  real  t£st  as  to  the  usefulness  and  speed  of  this  particular 
data  retrieval  system  since  our  data  base  will  vastly  exceed  any  data  base  previously  used.  When 
the  seismic  data  are  finally  entered  into  the  system,  the  user  will  be  able  to  search  the  hypo- 
centcr  file  or  the  tape  file  on  any  specific  parameter  or  set  of  parameters  and  conditions.  For 
example,  a  user  may  wish  a  list  of  all  tapes  that  contain  earthquakes  of  magnitude  5.0  or  greater, 
with  depths  less  than  60  km  from  the  Kuril  Islands  that  were  recorded  on  both  LASA  and  NORSAR 
tapes  in  our  library.  It  is  expected  that  the  LISTAR  system  will  produce  this  list  within  10  or  15 
minutes.  Such  a  manual  search  through  all  the  tape  listings  and  epicenter  listings  could  take  sev¬ 
eral  days. 

Construction  of  the  epicenter  data  file  is  now  under  way.  Initially,  the  epicenter  data  file 
will  contain  upwards  of  100,000  entries  obtained  from  as  many  reporting  agencies  as  will  be  con¬ 
venient.  It  is  anticipated  that  the  annual  increase  in  this  file  will  be  about  20,000  entries.  A  sec¬ 
ond  data  file,  our  tape  library  file,  will  initially  contain  about  8000  entries  and  be  increased  by 
1  500  entries  annually.  R  M  Shcppard 


t"A  User's  Guide  to  LISTAR,”  Lincoln  Manual  94,  Lincoln  Laboratory,  M.I.T.  (October  1970), 
DDC  AD-714108. 
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VII.  THE  INNER  CORE 


A.  DETECTION  OF  THE  PHASE  PKJKP 

Since  its  discovery  in  1935,  the  Earth’s  inner  core  has  been  assumed  to  be  a  solid,  but 
attempts  to  detect  seismic  waves  which  travel  through  it  as  shear  waves  have  always  failed. 

There  have  been  two  reasons  for  this  failure:  uncertainty  of  the  shear  velocity  in  the  inner  core 
and  thus  of  where  to  look  for  the  seismic  waves,  and  the  likelihood  that  any  such  waves  would  be 
extremely  small  and  might  be  masked  by  other,  larger  phases.  Solutions  to  both  problems  are 
now  available,  however,  so  a  search  for  the  phase  PKJKP  was  undertaken. 

In  the  first  place,  available  free  oscillation  data,  particularly  for  overtones  of  the  Earth's 
spheroidal  modes,  are  now  able  to  detect  the  inner  core’s  rigidity  and  to  make  a  rough  estimate 
of  its  value.  Thus,  it  is  possible  to  calculate  at  which  distances  and  travel  times  PKJKP  should 
be  largest.  Secondly,  large  seismic  arrays  such  as  LASA  and  NORSAR  can  be  used  to  separate 
PKJKP  from  "noise”  consisting  of  phases  with  different  apparent  velocities.  Normal  micro- 
scismic  noise  is  of  no  concern  here,  because  of  the  size  (m^  6.0)  of  the  events  used. 

Nine  large  earthquakes  of  various  depths  in  the  distance  range  from  70°  to  125°  (out  of  38 
events  examined  visually)  have  been  subjected  to  velocity  spectral  analysis  (Vespa,  see  Davies, 
et  al.^)  using  data  from  IASA.  Of  these,  five  had  phases  whose  apparent  velocities  could  be 
determined  reliably.  The  principal  source  of  confusion  in  this  respect  comes  from  sidelobes 
in  the  wavenumber  response  of  the  array,  so  care  has  been  taken  to  reject  signals  (some  prob¬ 
ably  real)  whenever  the  possibility  of  such  confusion  arises. 

Figure  VI I  - 1  shows  the  observed  travel  times  and  dT/dA  values  for  these  phases,  as  well 
as  two  theoretical  curves,  which  suggest  that  the  inner  core  has  a  shear  velocity  of  about 
2.95  ±0.1  (which  may  decrease  with  depth).  A  comparison  of  the  waveforms  of  PKKP  and  PKJKP 
consistently  indicates  that  PKJKP  is  depleted  in  frequencies  above  about  0.5  Hz,  but  that  other¬ 
wise  the  signal  shapes  are  remarkably  similar.  We  have  obtained  a  rough  estimate  of  the  Q  for 
the  J  phase  of  500  to  1000  by  comparing  the  waveform  of  PKJKP  with  that  of  PKKP  after  sub- 

lecting  it  to  an  attenuation  filter.^  ~  T  i 

J  6  B.  R.  Julian 

D.  Davies 

R.  M.  Sheppard 


B.  STRUCTURE  OF  THE  INNER-CORE  BOUNDARY 

A  very  interesting  prospect  arises  from  the  work  reported  in  Sec.  I  by  Davies  —  a  means  of 
determining  the  fine  structure  of  the  inner-core  boundary.  If  we  take  PcP  as  a  source  function 
and  deconvolve  it  from  PKiKP,  then  the  resulting  time  function  should  be  a  measure  of  the  re¬ 
flectivity  at  the  inner-core  boundary  as  a  function  of  depth.  The  assumption  is  that  the  outer 
core  has  a  very  sharp  boundary  —  which  seems  justifiable  from  the  strong  similarity  between  PcP 
and  PKKKP.  Figure  Vll-2  shows  the  transfer  function  for  inner-core  reflection  as  a  function 


t  D.  Davies,  E.  J.  Kelly  and  J.  R.  Filson,  "Vespa  Process  for  the  Analysis  of  Seismic  Signals," 
Nature  232,  8-13  (1971). 

%  E.  W.  Carpenter,  "Absorption  of  Elastic  Waves  —  an  Operator  for  a  Constant  Q  Mechanism," 
Atomic  Weapons  Research  Establishment  Report  0-43/66,  Her  Majesty's  Stationery  Office, 
London  (1  966). 
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Section  VI! 


of  time.  There  is  obvious  structure  over  an  interval  of  at  least  2  sec,  from  which  we  infer  that 
the  inner-core  boundary  is  not  sharp  to  better  than  10km.  This  is  an  early  conclusion  and  sub¬ 
ject  to  change  when  more  data  have  been  examined.  Visual  inspection  of  the  signal  following 
PKiKP  suggests  that  there  may  be  structure  for  a  further  100  km  into  the  inner  core. 

D.  Davies 
C.  W.  Frasier 


MO  ITO  250  230 

DELTA  (  d*g) 


Fig.  VI 1-1.  Observed  travel  times  and  dT/dA  values  for  PKJKP,  with  theoretical 
curves  for  inner-care  velocities  (assumed  constant)  of  2. 9  and  3.0  km/sec.  Data 
are  alsa  shown  far  various  branches  of  phase  PKKP. 
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Fig.  VI 1-2.  A  measure  of  reflectivity  as  a  function  of  time  at  inner-core  boundary. 
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